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FOREWORD 
The  Analytical  Chemistry  Division  was  established  as  a 
separate  division  at  the  National  Bureau  of  Standards  on 
September  1,  1963,  and  became  part  of  the  Institute  for 
Materials  Research  in  the  February  1,  196H  reorganization. 
It  consists  at  present  of  nine  Sections  and  about  100  tech- 
nical personnel  encompassing  some  45  different  analytical 
competences  from  activation  analysis  and  atomic  absorption 
to  vacuum  fusion  and  x-ray  spectroscopy.   These  competences, 
and  in  turn  the  Sections  which  they  comprise,  are  charged 
with  research  at  the  forefront  of  analysis  as  well  as  aware- 
ness of  the  practical  sample,  be  it  standard  reference 
material  or  service  analysis.   In  addition  it  is  their 
responsibility  to  inform  others  of  their  efforts. 

Formal  publication  in  scientific  periodicals  is  a  highly 
important  output  of  our  laboratories.   In  addition,  however, 
it  has  been  our  experience  that  informal,  annual  summaries  of 
progress  describing  efforts  of  the  past  year  can  be  very  valu- 
able in  disseminating  information  about  our  programs.   A  word 
is  perhaps  in  order  about  the  philosophy  of  these  yearly 
progress  reports.   In  any  research  program  a  large  amount  of 
information  is  obtained  and  techniques  developed  which  never 
find  their  way  into  the  literature.   This  includes  the 
"negative  results"  which  are  so  disappointing  and  unspectacu- 
lar but  which  can  often  save  others  considerable  work.   Of 
inportance  also  are  the  numerous  small  items  which  are  often 
explored  in  a  few  days  and  which  are  not  important  enough  to 
warrent  publication  -  yet  can  be  of  great  interest  and  use 
to  specialists  in  a  given  area.   Finally  there  are  the  exper- 
imental techniques  and  procedures,  the  designs  and  modif ications 
of  equipment,  etc.,  which  often  require  months  to  perfect  and 
yet  all  too  often  must  be  covered  in  only  a  line  or  two  of  a 
journal  article. 
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Thus  our  progress  reports  endeavor  to  present  this  in- 
formation which  we  have  struggled  to  obtain  and  which  we 
feel  might  be  of  some  help  to  others.   Certain  areas  which 
it  appears  will  not  be  treated  fully  in  regular  publications 
are  considered  in  some  detail  here.   Other  results  which  are 
being  written  up  for  publication  in  the  journal  literature 
are  covered  in  a  much  more  abbreviated  form. 

At  the  National  Bureau  of  Standards  publications  such 
as  these  fit  logically  into  the  category  of  a  Technical  Note. 
We  plan  to  issue  these  summaries  for  all  of  our 
Sections.   The  following  is  the  sixth  annual  report  on 
progress  of  the  Radiochemical  Analysis  Section. 

W.  Wayne  Meinke,  Chief 
Analytical  Chemistry  Division 
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ABSTRACT 

This  is  the  sixth  summary  of  progress  of  the  Radiochem- 
ical Analysis  Section  of  the  Analytical  Chemistry  Division 
at  the  National  Bureau  of  Standards . 

The  section's  effort  comprises  five  major  areas: 
Mossbauer  spectroscopy,  nuclear  chemistry,  nuclear  instru- 
mentation, radioisotope  tracer  techniques  and  the  application 
of  statistics  in  nuclear  and  analytical  chemistry. 

Low  temperature  devices  for  Mossbauer  spectroscopy  are 
a  subject  of  continuing  Interest,  and  this  year's  effort  is 
centered  around  the  design  of  variable  temperature  devices. 
A  high  temperature  device  is  also  described. 

Preliminary  data  on  the  detection  of  conversion  elec- 
trons show  the  potential  for  Mossbauer  spectroscopy  in  the 
field  of  chemical  bonding  in  surfaces. 

Further  exploration  of  the  potential  of  Mossbauer 
spectroscopy  for  the  measure  of  structural  parameters  of 
nickel  compounds  have  been  made . 

Preliminary  data  are  presented  on  the  structure  of  iron 
doped  rutile  (TiOp)  which  is  of  significance  to  the  deter- 
mination of  mineral  structures  found  on  the  moon.   Theoret- 
ical predictions  of  Mossbauer  spectra  resulting  from  electric 
field  gradients  due  to  crystal  field  interactions  are  pre- 
sented.  Evaluation  of  the  basic  integrals  using  the  self- 
consistent  Hartree-Fock  method  is  given  for  ions  of  the  first 
row  of  transition  elements. 
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Precautions  to  be  taken  in  the  analysis  of  gamma-ray 
spectra  are  presented  when  model  errors  exist  and  when  a 
gain  shift  is  detected. 

A  system  for  gas  analysis  which  utilizes  gas  chroma- 
tography is  described.  Its  application  to  nuclear  cross 
section  studies  is  also  presented.  A  method  for  setting 
limits  on  sample  heterogeneity  using  a  specific  analysis 
method  provides  a  useful  indicator  for  evaluation  of  the 
adequacy  of  standard  reference  materials . 

An  interesting  presentation  of  detection  and  measure- 
ment in  the  femtogram  region  is  given. 

The  complete  circuitry  for  the  optical  interfero- 
metric  Mossbauer  spectrometer  is  presented  in  this  report. 

Quantitative  analysis  for  molybdenum  and  iron  using 
radioisotope  dilution  to  a  precision  that  rivals  most 
other  techniques  is  described. 

James  R.  DeVoe ,  Chief 
Radiochemical  Analysis  Section 
Analytical  Chemistry  Division 

Key  Words : 
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Mossbauer  spectrometry;  nickel-6l;  rutile;  structure 
analysis;  counting  statistics;  heterogeneity;  proportional 
detector;  conversion  electrons;  electric  field  gradient 
tensor;  internal  magnetic  field;  self-consistent  crystal 
field;  calculations;  interferometric ;  Mossbauer  spectrom- 
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1.   INTRODUCTION 

The  Radiochemical  Analysis  Section  develops  measurement 
techniques  for  qualitative,  quantitative  and  structural  anal- 
ysis of  materials  through  the  use  of  radioisotopes.   These 
techniques  require  a  very  interesting  and  stimulating  blend 
of  knowledge  ranging  from  nuclear  and  atomic  physics  through 
nuclear  chemistry  and  physical  chemistry  to  descriptive 
chemistry . 

It  Is  essential  that  the  methods  developed  be  applica- 
ble in  the  sense  that  they  can  be  used  successfully  for 
materials  of  interest  to  science  and  industry.   Therefore, 
these  methods  are  applied  to  the  analysis  of  NBS  Standard 
Reference  Materials  where  considerable  cross  checking  of 
analytical  techniques  is  required  throughout  the  process 
of  certification. 

In  order  to  define  more  accurately  the  scope  of  our 
work,  abstracts  of  publications  appear  in  this  report. 
Brief  descriptions  of  studies  nearing  completion  are  in- 
cluded when  it  is  our  intent  to  submit  this  work  elsewhere. 

A  roster  of  the  groups  in  this  section  is  listed  in 
part  5.   The  National  Bureau  of  Standards  has  several  pro- 
grams whereby  a  scientist  from  the  United  States  or  abroad 
may  work  in  our  laboratories  for  one  or  two  years .   It  is 
hoped  that  by  utilizing  these  programs  the  section  will  be 
able  to  perpetuate  a  stimulating  environment. 

In  order  to  specify  adequately  the  procedures,  it  has 
been  necessary  occasionally  to  identify  commercial  materials 
and  equipment  in  this  report.   In  no  case  does  such  identi- 
fication imply  recommendation  or  endorsement  by  the  National 
Bureau  of  Standards,  nor  does  it  imply  that  the  material  or 
equipment  identified  is  necessarily  the  best  available  for 
the  purpose. 

James  R.  DeVoe ,  Chief 
Radiochemical  Analysis  Section 
Analytical  Chemistry  Division 
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2.   MOSSBAUER  SPECTROSCOPY 

A .  Introduction 

Our  efforts  remain  centered  about  improving  instrumen- 
tation, spectral  interpretation,  and  of  most  importance,  the 
establishment  of  suitable  standard  materials  for  the  purpose 
of  improving  the  accuracy  of  spectral  data.   The  basic  de- 
sign for  the  constant  acceleration  interferometer-controlled 
Mossbauer  spectrometer  was  completed.   Refinements  made  on 
the  system  have  resulted  in  a  spectrometer  with  exceptionally 
fine  qualities.   A  resonant  detector  was  made  operational. 
Its  performance  has  provided  much  information  on  the  para- 
meters of  design  that  must  be  considered  to  optimize  the 
overall  efficiency  of  such  detectors.   Another  very  signifi- 
cant enhancement  of  our  capability  is  the  availability  of 
several  versatile  systems  which  allow  measurement  of  Mossbauer 
spectra  at  temperatures  ranging  from  4  to  1000  K. 

Another  major  accomplishment  relates  to  improvements 
made  on  the  nonlinear  least  squares  fitting  procedures 
which  incorporate  a  modified  convergence  criterion  and  the 
use  of  basic  structural  parameters.   This  work,  coupled  with 
mathematical  procedures  for  calculating  electric  fields  from 
single  crystal  data,  has  enhanced  our  capability  for  inter- 
pretation of  spectra. 

Therefore,  the  group  is  now  in  the  excellent  position 
of  having  done  its  "homework",  and  it  now  enters  a  period  of 
anticipated  high  productivity.   Heading  the  list  are  standard 
materials  in  the  form  of  iron  metal  for  velocity  calibration 
and  bench  marking  isomer  shift,  a  recalibration  of  sodium 
nitroprusside  for  the  isomer  shift  standard,  and  a  calibra- 
tion of  dimethyl  tin  difluoride  for  the  isomer  shift  standard. 

B.  Instrumentation 

1 .   Helium  Temperature  Cryostat  for   Ni  Experiments 
The  low  recoil-free  fraction  (f)  observed  in  nickel 
compounds  requires  Mossbauer  measurements  to  be  made  at 
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liquid  helium  temperatures.   Since  the  attenuation  of  67.^ 
keV  gamma  radiation  by  borosilicate  glass  is  small,  a  simple 
helium  cryostat  was  designed  for  cooled  source  and  absorber, 
with  rapid  sample  changing  capability.   The  cryostat  is 
shown  in  Figure  1.   The  glass  Dewar  (A)  is  suspended  from  an 
aluminum  support  plate  (B)  by  a  split  clamping  ring  (C)  and 
Teflon*  gasket  (D) .   The  Mossbauer  drive  is  mounted  on  the 
aluminum  support  plate  with  an  O-ring  seal  (E),  and  the  source- 
absorber  holders  (P)  are  rigidly  attached  to  the  drive  (G), 
such  that  the  drive  and  source-absorber  assembly  can  be  lifted 
out  of  the  cryostat.   The  source  and  absorber  assembly  is 
shown  in  Figure  1-a.   The  1/8  inch  diameter  tubing  is  rigid 
enough  to  prevent  vibration  and  has  low  thermal  conduction. 
The  source  is  stationary  and  the  absorber  is  mounted  on  a 
flexure  plate  (I),  and  coupled  to  the  drive  by  a  yoke  (J) 
and  1/8  inch  stainless  steel  tubing  (K).   Source  or  absorber 
can  be  changed  in  5  minutes.   The  67.^  keV  gamma-ray  atten- 
uation is  15%,    and  liquid  helium  consumption  is  10  h/liter. 

(D.  Sweger  and  J.  J.  Spijkerman) 

2 .   Simple  Variable  Temperature  Cryostat  for 
Temperature  Range  of  100  to  373  K 

The  cryostat  described  in  Reference  [1],  p25,  has  been 
modified  for  variable  temperature  use,  by  inserting  a  carbon 
cloth  heater  between  the  heat  sink  and  the  sample  mounting 
block.   The  1  inch  by  1.5  inch  carbon  cloth  is  clamped  between 
two  copper  foil  electrodes,  and  insulated  by  1/16  inch  Teflon* 
plates,  (see  Figure  2).   The  heater  is  operated  by  a  6.3  volt 
filament  transformer  connected  to  a  temperature  controller  or 
variable  transformer  for  manual  control.   For  manual  control, 
the  temperature  is  stabilized  after  30  minutes  to  within  ± 
2°C  over  a  period  of  6  hours  and  is  measured  by  a  copper 
constantan  thermocouple. 

(P.  R.  Stampfl  and  J.  J.  Spijkerman) 


*For  disclaimer  of  equipment  and  materials  see  last  paragraph 
of  the  Introduction. 


Figure  la 


Figure  1.   Low  temperature  cryostat.   Fig.  la.  Source 

absorber  assembly  for  low  temperature  cryostat 
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Figure    2.      Heater   for   source-absorber   assembly    of    simple 
helium  cryostat. 
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3 .  Furnace  for  Absorber  Temperature  Studies  from  300 
to  1000  K 

Conventional  type  furnaces  are  not  very  well  suited  for 
Mossbauer  work  due  to  their  size  and  difficulty  in  sample 
mounting.   By  reducing  the  length  of  the  furnace,  the  time 
to  accumulate  data  can  be  greatly  reduced   thereby  also 
reducing  the  requirement  for  long-time  temperature  stability. 
However,  the  temperature  uniformity  across  the  sample  is 
affected.   This  can  be  eliminated  if  the  sample  is  placed  in 
contact  with  the  heater,  as  shown  in  Figure  3-   The  sample 
is  clamped  between  two  BeO  disks,  2  mm  thick,  which  are  in 
direct  contact  with  two  carbon  cloth  heaters.   Thin  sheets 
of  nickel  plated  mica  provide  a  radiation  shield,  and  alu- 
mina rings  are  used  for  support  and  temperature  insulation. 
Including  the  vacuum  jacket,  the  furnace  has  a  2  inch  outer 
diameter  and  is  1.5  inches  long.   Temperatures  of  1000  K 
are  obtained  and  were  measured  with  a  chromel-alumel 
thermocouple . 

(J.  J.  Spijkerman) 

4 .  Conversion  Electron  Proportional  Counter  for 
Spectral  Measurement  with  Backscattering  Geometry 

Experiments  with  the  pancake  detector  described  in 
Reference  [1],  p.  28,  for  detecting  conversion  electrons 
showed  that  spectra  with  16%  effect  above  baseline  could 
be  obtained.   A  new  detector  was  designed  and  optimized  for 
conversion  electron  detection  (see  Figures  4  and  5).   The 
detector  has  two  0.001  inch  wire  electrodes,  spaced  1/4  inch 
apart.   The  detector  operates  at  1200  to  1300  V  and  uses  10$ 
CH^-90%  He  flow  gas.   The  effective  penetration  depth  into 
the  sample  for  Mossbauer  analysis  was  determined  from  the 
spectra  (shown  in  Figures  6  and  7)  of  iron  evaporated  on 
stainless  steel  foils.   The  contribution  from  a  layer  at  a 
distance  (y)  from  the  surface  to  the  Mossbauer  backscattering 
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Figure  3.   Exploded  diagram  of  high  temperature  furnace 
for  absorber. 
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3 •   Furnace  for  Absorber  Temperature  Studies  from  300 
to  1000  K 

Conventional  type  furnaces  are  not  very  well  suited  for 
Mossbauer  work  due  to  their  size  and  difficulty  in  sample 
mounting.   By  reducing  the  length  of  the  furnace,  the  time 
to  accumulate  data  can  be  greatly  reduced   thereby  also 
reducing  the  requirement  for  long-time  temperature  stability. 
However,  the  temperature  uniformity  across  the  sample  is 
affected.   This  can  be  eliminated  if  the  sample  is  placed  in 
contact  with  the  heater,  as  shown  in  Figure  3.   The  sample 
is  clamped  between  two  BeO  disks,  2  mm  thick,  which  are  in 
direct  contact  with  two  carbon  cloth  heaters.   Thin  sheets 
of  nickel  plated  mica  provide  a  radiation  shield,  and  alu- 
mina rings  are  used  for  support  and  temperature  insulation. 
Including  the  vacuum  jacket,  the  furnace  has  a  2  inch  outer 
diameter  and  is  1.5  inches  long.   Temperatures  of  1000  K 
are  obtained  and  were  measured  with  a  chromel-alumel 
thermocouple . 

(J.  J.  Spijkerman) 
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showed  that  spectra  with  16%  effect  above  baseline  could 
be  obtained.   A  new  detector  was  designed  and  optimized  for 
conversion  electron  detection  (see  Figures  4  and  5).   The 
detector  has  two  0.001  inch  wire  electrodes,  spaced  1/4  inch 
apart.   The  detector  operates  at  1200  to  1300  V  and  uses  10% 
CHn-90%  He  flow  gas.   The  effective  penetration  depth  into 
the  sample  for  Mossbauer  analysis  was  determined  from  the 
spectra  (shown  in  Figures  6  and  7)  of  iron  evaporated  on 
stainless  steel  foils.   The  contribution  from  a  layer  at  a 
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intensity  for  resonant  velocity  can  be  represented  by 


arc    cos    y/R 


/°°           _       „                               -,  rare    cos    y 

dEe   naJ,a    °lHy    ' -^ — ^ J  sin    6    d    9 


(1) 


E   +  Ts  /4  0 


Where  9  =  angle  of  scattering  from  the  normal  to  the  surface 
C  =  constant 
y  =  absorption  coefficient  for  the  incident  14.4  keV 

gamma  ray 
R  =  the  range  of  an  8  keV  electron  in  iron. 
E  =  Mossbauer  gamma-ray  energy 
n  =  number  of  resonant  atoms  per  cm   in  the  absorber. 

f   =  fraction  of  recoilless  absorption  in  absorber. 

a  v 

0-,  u   =  cross  section  for  resonant  absorption. 
T  =  half  width  of  source. 
Since  R  is  very  small,  the  attenuation  of  the  14.4  keV  gamma 
radiation  can  be  neglected,  and  y<R, 

T(y)  =C'|   (l  -  |)  (2) 

and  the  backscattering  intensity  I(y)  is 

I(y)  =   T(Y)dy  =  C  (  2L.  _  L  (3) 


\2R   3R2  / 


normalizing,  such  that  I(R)  =  1 


(4) 


The  backscattering  amplitude  I(y)  as  a  function  of  depth  into 
the  sample  is  shown  in  Figure  o,  for  R=1200  A.   The  Mossbauer 
backscatter  intensity  for  600  A  iron  on  stainless  steel  foil 
(Figure  7)  was  65% .   From  Figure  8,  this  would  represent  the 

o 

effective  penetration  depth  of  approximately  1000  A.   However, 

2950  A  of  iron  (Figure  6)  on  stainless  steel  foil  gave  96%, 
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Figure  8.   Backscatter  amplitude  vs.  depth  of  penetration 
into  sample. 

with  h%   of  the  intensity  originating  from  the  stainless  steel 
foil.   This  indicates  that  there  is  a  serious  difficulty  with 
the  model  for  the  range  of  8  keV  electrons  in  iron.   Further 
experimental  work  is  in  progress. 

(J.  J.  Spijkerman) 

5 .   Resonant  Detector 

a .   Detection  of  Mossbauer  Events  by  Internal 

Conversion  Electrons:   Preliminary  Theoretical 
Studies 
The  resonance  detection  should  provide  unambiguous 
recognition  of  a  Mossbauer  event  characterized  by  the  re- 
emission  of  the  gamma  photon  one  time  out  of  ten,  or  the 
emission  of  an  internal  conversion  electron  and  its  attendant 
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intensity  for  resonant  velocity  can  be  represented  by 


/°°                    -.        /-arc  cos  y/R 
dEe  na*a  °lhy    '-= — ^-= /  sin  9  d  9 

1  E^  +  T^  /4  J7 
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Where  9  =  angle  of  scattering  from  the  normal  to  the  surface 

C  =  constant 

u  =  absorption  coefficient  for  the  incident  lH.k   keV 

gamma  ray 

R  =  the  range  of  an  8  keV  electron  in  iron. 

E  =  Mossbauer  gamma-ray  energy 
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radiation  can  be  neglected,  and  y<R, 

T(y)  =  0'  |  (l  -  |)  (2) 

and  the  backscattering  intensity  I(y)  is 

I(y)  =   T(Y)dy  =  C  (  Z_  _  £_  }  (3) 


\  2R   3R2  / 
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The  backscattering  amplitude  I(y)  as  a  function  of  depth  into 
the  sample  is  shown  in  Figure  0,  for  R=1200  A.   The  Mossbauer 
backscatter  intensity  for  600  A  iron  on  stainless  steel  foil 
(Figure  7)  was  65%.   From  Figure  8,  this  would  represent  the 
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foil.   This  indicates  that  there  is  a  serious  difficulty  with 
the  model  for  the  range  of  8  keV  electrons  in  iron.   Further 
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5 .   Resonant  Detector 

a .   Detection  of  Mossbauer  Events  by  Internal 

Conversion  Electrons:   Preliminary  Theoretical 
Studies 
The  resonance  detection  should  provide  unambiguous 
recognition  of  a  Mossbauer  event  characterized  by  the  re- 
emission  of  the  gamma  photon  one  time  out  of  ten,  or  the 
emission  of  an  internal  conversion  electron  and  its  attendant 
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x-ray  nine  times  out  of  ten.   Of  these  three  indicators  of 
Mossbauer  events,  the  photons  may  be  confused  with  similar 
photons  entering  the  detector  from  the  source  and  absorber, 
and  thus  are  hardly  unambiguous.   The  internal  conversion 
electron,  on  the  other  hand,  cannot  be  confused  with  electrons 
originating  outside  of  the  counter,  since  these  would  never 
reach  the  inside  of  the  counter,  but  they  can  be  confused 
with  photoelectrons  and  Compton  electrons  created  in  the 
Mossbauer  gamma  ray  and  the  higher  energy  gamma  radiation 
from  the  source. 

The  following  paragraphs  and  charts  outline  some  very 
preliminary  theoretical  studies  designed  to  predict  the 
contributions  of  the  various  electron  producing  processes . 
For  this  first  study,  several  simplifying  assumptions  and 
approximations  have  been  made : 

1.  Only  internal  conversion  (Mossbauer)  and  photo- 
electrons  have  been  considered.   This  is  a  good  approximation 
of  the  lU . 4  keV  gamma  ray,  but  not  for  those  of  122  and  136 
keV  which  should  be  Compton  scattered  to  roughly  the  same 
degree  as  photoabsorbed. 

2.  A  natural  abundance  (in  Fe-57)  thick  iron  foil 
has  been  assumed  as  the  in-detector  absorber,  and  the  effect 
of  the  detector  walls  has  been  neglected. 

3.  Auger  electrons  and  secondary  interaction  elec- 
trons have  been  neglected. 

4.  Discrete  gamma-ray  energies  have  been  assumed, 
ignoring  natural  line  broadening  and  Compton  smearing  of  the 
high  energy  gammas  in  any  external  absorber  which  may  be 
present . 

5.  Only  electrons  traveling  the  shortest  distance 
from  the  point  of  origin  to  the  iron  surface  have  been  con- 
sidered.  Angular  distributions  will  be  included  in  future 
calculations . 
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6.   A  "most  probable"  energy  loss,  as  opposed  to  a 
distribution  of  energy  losses,  has  been  assumed  for  each 
electron  as  a  result  of  deceleration  in  the  foil. 

The  number  of  electrons  created  over  a  small  increment 
of  thickness  at  depth  'x1  into  the  material  is  given  by 

N   (x)  =  a  I  exp   -(a   +  a.jN^-,  x  (c) 

p        poKp     M   F  v^ 

for  photoelectrons  and 


NM(x)  =  aMIoexp   "(ap  +  aM)NF  x'  a/(]  : 


] 


(6) 


for  internal  conversion  electrons.   Here  I   is  the  relative 

o 

intensity  of  the  parent  gamma,  a   and  aM  are  the  photo- 
electric and  Mossbauer  cross  section,  Nn  is  the  number  of 

3  F 

atoms  per  cmJ  in  the  foil,  and  a  is  the  internal  conversion 

coefficient.   These  equations  are  used  independently  for  each 

of  the  incident  gammas,  with  the  appropriate  intensities  and 

cross  sections. 

The  most  probable  energy  loss  for  an  electron  of 

kinetic  energy  E,  escaping  from  a  depth  of  'x'  is  taken  to  be 

(4  E   £y 
log  -=-^ — =-  -  g2  +  0.37  )  (7) 

i  (i-e  ) 


32   =    1   -    hll   keV/(5H   keV  +   E    )  1 


K    =   5.613    2   x   10    10   keV   cm"1 

where  I  is  the  average  ionization  potential  for  the  foil. 

The  resultant  energy  carried  by  the  electron  escaping  from  the 

foil  is  simply  E   -  AT  ,  where  E.  =  E   -  B,  and  B  is  the 

k.      p  k     y 

electron  binding  energy. 

Table  I  and  Figure  9  show  the  calculated  resultant 
energy  and  relative  intensities  of  Mossbauer  and  photoelec- 
trons emerging  from  a  specified  depth  within  an  iron  foil  for 
the  14.4,  122,  and  136  keV  gamma  rays  and  cobalt-57. 

(J.  C.  Travis) 
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Table  1.   Photo-electron  and  conversion  electron  production 
in  iron  metal . 

Relative  Intensities 


Depth 

Photo- 

Mossbauer 

Resultant 

Gamma 

(A) 

electrons 

conversion 

energy 

energy 

electrons 

(keV) 

(keV) 

500 

54510.9 

498459 

6.57612 

14.4 

1500 

51909.6 

474673. 

5.38439 

14.4 

2500 

49432.5 

452021. 

4.04075 

14.4 

3500 

47073.6 

430451. 

2.61123 

14.4 

4500 

44827.3 

409910. 

1.12123 

14.4 

25000 

83582.7 

0 

111.097 

122 

75000 

80628.7 

0 

102.703 

122 

125000 

77779.1 

0 

93.6968 

122 

175000 

75030.3 

0 

84.3441 

122 

225000 

72378.6 

0 

74.7475 

122 

275000 

69820.6 

0 

64.9623 

122 

325000 

67353. 

0 

55.0232 

122 

375000 

64972.6 

0 

44.954 

122 

425000 

62676.3 

0 

34.7723 

122 

475000 

60461.2 

0 

24.4913 

122 

525000 

58324.4 

0 

14.1215 

122 

575000 

56263.1 

0 

3.67144 

122 

25000 

7218.94 

0 

125.358 

136 

75000 

7044.06 

0 

117.598 

136 

125000 

6873.^2 

0 

109.274 

136 

175000 

6706.91 

0 

100.632 

136 

225000 

6544.43 

0 

91.7652 

136 

275000 

6385.9 

0 

82.7248 

136 

325000 

6231.2 

0 

73.5428 

136 

375000 

6080.25 

0 

64.2411 

136 

425000 

5932.95 

0 

54.8359 

136 

475000 

5789.23 

0 

45.3393 

136 

525000 

5648.98 

0 

35.761 

136 

575000 

5512.14 

0 

26.1088 

136 

625000 

5378.61 

0 

16.3892 

136 

675000 

5248.31 

0 

6.60758 

136 
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Figure  9.   Energies  and  relative  numbers  of  conversion  and 
photoelectrons  vs  depth  into  iron  foil. 


b .   Design  of  Resonant  Detector 

The  theoretical  study  shows  that  the  detection  of 
conversion  electrons  with  optimum  electron  collection  in 
the  detector  can  give  a  maximum  signal/noise  ratio  of  500, 
using  100%  enriched  iron  stainless  steel.   The  conversion 
electron  counter  described  on  page  6  with  a  single  50% 
enriched  0.001  inch  stainless  steel  foil  was  used  to  deter- 
mine the  signal/noise  ratio.   The  spectrum  shown  in  Figure 
10  gave  a  S/N  ratio  of  3,  using  only  lower  level  energy 
discrimination.   Spectra  were  taken  for  several  single 
channel  settings  to  determine  the  electron  energy  distri- 
bution in  the  counter  for  the  8  keV  conversion  electrons . 
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Figure  11 


3        4 
E  (volts) 


Signal  to  noise  (S/N)  ratio  as  a  function  of 
single  channel  setting  using  conversion  electron 
detection  from  stainless  steel. 


The  results  shown  in  Figure  11,  indicate  that  the  energy 
distribution  of  the  electrons  escaping  from  the  foil  is 

shifted  toward  lower  energies.   The  signal/noise  can  be 
improved  by  using  a  multilayer  detector,  with  1000  A  stain- 
less steel  foils  containing  highly  enriched  iron-57.   This 
will  also  increase  the  counting  efficiency,  and  lower  the 
photoelectron  noise  originating  from  the  122  keV  precursor 
gamma  ray.   Please  see  acknowledgements. 

(J.  J.  Spijkerman  and  J.  C.  Travis) 
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C .   SRM  for  the  Isomer  Shift  of  Tin  Compounds 

After  a  long  search  for  a  suitable  tin  Mossbauer 

isomer  shift  standard,  a  stable  compound  with  a  relatively 

large  quadrupole  splitting  was  obtained.   Dimethyl  tin 

difluoride  has  a  quadrupole  splitting  of  about  4.21  mm/sec 

and  a  narrow  line  width.   Calibration  of  this  material  will 

be  completed  in  November  1969 .   The  spectrum  for  a  100 

mg/cm   sample  is  shown  in  Figure  12.   Several  plastic 

encapsulated  samples  of  various  thicknesses  supplied  by  Dr. 

Beasley  of  Austin  Science  Associates  were  analyzed  with  the 

results  shown  in  Table  2.   The  optimum  sample  thickness 

2 
would  be  100  mg/cm  . 

Table  2.   Percent  effect  vs.  thickness  of  dimethyl  tin 
difluoride . 

Sample  Line  width 

Thickness 

z   2 

mg/cm 

50 
100 
150 


Percent 

effect 

mm/sec 

FWHM 

Peak  1 

Peak  2 

Peak  1 

Peak  2 

4.2 

7.5 

10.6 

3.5 
6.3 
9.4 

.87 
.90 
.96 

.87 
.89 
.94 

D .   Constrained  Fitting  of  Complex  Mossbauer  Spectra 

The  concepts  described  earlier  [1],  for  the  fitting 
of  unresolved  nickel-6l  Mossbauer  spectra  using  the  OMNITAB 
programming  language  have  been  adapted  to  a  nonlinear  least 
squares,  FORTRAN  curve  fitting  program.   In  addition,  the 
capability  of  fitting  multisite  iron-57  spectra  has  been 
added.   The  new  program  which  uses  the  algorithm  devised  by 
D.  W.  Marquart  [2]  provided  a  balance  between  the  Gauss- 
Newton  method,  used  in  the  TN451  [1]  and  TN404  [3]  (PARLORS) 
programs,  and  the  method  of  steepest  descent.   Inclusion 
of  the  steepest  descent  method  stabilizes  the  calculation  in 
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regions  of  the  parameter  space  where  the  first  order  Taylor 
series  expansion  of  the  Gauss-Newton  method  is  not  accurate. 
Although  nonconvergence  is  thus  eliminated,  caution  is  still 
required  in  making  initial  estimates  to  avoid  convergence 
to  a  false  minimum. 

The  program  is  currently  undergoing  revision  to 
improve  the  input  and  output  and  to  decrease  the  running 
time.   Program  decks  and  documentation  will  be  made  available 
'!  at  a  later  date  . 

The  figures  included  here  illustrate  the  use  of  the 
program  in  resolving  complex  spectra.  Figures  13  and  14  are 
twelve  line  magnetic  nickel-6l  spectra  with  and  without  quad- 
rupole  perturbation,  respectively.  Magnetic  iron-57  spectra 
of  systems  with  one  and  three  sites,  respectively,  are  shown 
in  Figures  15  and  16.  Figure  17  shows  iron-57  with  two  mag- 
netic sites  and  one  single  line  site,  and  Figure  18  shows 
a  nonmagnetic  iron  spectrum  fitted  with  two  quadrupole  split 
sites  and  an  unsplit  site. 

The  number  of  sites  to  be  fitted  is  governed  by  the 
number  of  input  estimate  cards.   Each  estimate  card  (one 
for  each  site)  contains  four  input  parameters;  1)  centroid 
of  the  pattern,  2)  halfwidth,  3)  magnetic  splitting  separa- 
tion (channels)  of  the  outermost  peaks  ,  and  4)  quadrupole 
splitting.   The  nature  of  each  site  (magnetic,  quadrupole 
split,  etc.)  is  indicated  by  the  input  estimates.   For 
example,  for  iron-57,  a  singlet  is  fitted  if  the  initial 
magnetic  field  and  quadrupole  splitting  estimates  are  zero, 
a  symmetric  doublet  is  fitted  if  only  the  magnetic  field 
estimate  is  zero,  and  a  six  line  pattern  is  fitted  other- 
wise.  The  six  line  pattern  is  constrained  to  have  relative 
areas  of  3:2:1:1:2:3,  with  the  half  widths  and  amplitudes 
modified  by  a  thickness  parameter.   Quadrupole  splitting  is 
included  to  first  order  for  a  nonzero  initial  estimate,  or 
ignored  for  an  initial  estimate  of  zero. 

21 


o 

ID 


(NJ 

O 
O 


II 
LU 


<X 
CD 


© 


o 

CO 

ro 


o 

(NJ 

in 

—) 

o 


ID 

a:    i 

i— 
<_j 

LU 
Q_ 
CO 


e© 


o 

o 

00 

. 

W 

CM 

• 

o 

^r 

■< 

■ 

n 

■P 

cd 

■=r 

O 

o 

C\J 

00 

"n 

■* 

O 

•H 

S 

Cm 

o 

O 

CM 

• 

(O 

3 

\ 

M 

E 

■P 

3E 

o 

0 

no 

a 

-^ 

ra 

I— 

M 

?H 

o 

^° 

o— I 

•LU 

CD 
ctf 

"?> 

CO 

:o 

S 

o 

H 

id 

CD 

V 

Ji4 

i 


0»*20I    09*101    08*001    00*001     02*66     0t*86 

A1ISN31NI    3AUHl^y 


09-Z.6 


08*96 


i 


PO 


CD 

M 

bO 
•H 


22 


o 
o 

^ 

r  © 

©    A 

in 

© 

a      ©« 

o 
o 

© 
e 

j>© 

© 
© 

% 

i       e 

o 

o 

r> 

©0 

®\©            © 

*«n 

CO 

oo 

o 
o 

(M 
II 

S|^>^ 

eg 

LU 

e       — .S 

z 

-»^^           © 

1—1 

£  ;® 

0    */        e 

o 

_l 

o 

LU 

■^ 

CO 

cr 

CD 

o 



o 
o 



IO 

\J}      © 

"o 

h- 

o 
o 

- 

© 

€T 

• 

1 

*"^"©© 

©               Qa^-^®"^ 

o 

G         ^CT® 

o 

m 

e>/®^ 

-CM 

LD 

*2^ 

•-H 

/?© 

^H 

*/               © 

o 

<n 

e 

y  G 

o 

• 

•  z.  ••• 

-en 

i 

O 

0 

©/        e 

z 

© 

e°° 

© 

i> 

f        G    © 

o 
o 

• 

CXL 

4       e     © 

.-* 

o 

©  © 

]©e  ° 

LU 

©f 

Q_ 

© 

o 

(O 

o 

08*00 

I 

OfOOT 

00 

•001 

09*66            02*66 

08*86            Of  86 

00*8 

6' 

CO 


o 
o 


"T> 


CM 


■P 

cti 

C\J 

CQ 
•H 

O 


3 
in 
-P 
O 
CD 
0- 
CO 

*H 

CD 

cd 

CO 
CO 

:0 


0) 

J* 
o 

■H 


0 

•H 


A1ISN31NI    3AIlB13y 


23 


o 
o 


CD 
■P 

•H 
CO 


o 

o 

o 

<o 

•H 

-p 

(D 

C 

hO 

o 

cd 

o 

£ 

• 

■* 

CD 

O 

-p 

•H 


3 

P 

o 

CD 

a 

CO 
CD 

CO 
CO 

:o 


o 

H 


tH 

CD 
U 

•H 


00 ' 00  I 


O0'B6 


00-96 


00't6 

JUISN31NI 


00 '26  00-1 

3AIl«13iJ 


oo- aa 


00-99 


24 


o 
o 


CD 
Ui 

CO 
00 

o 

CM 

II 


UJ 

to 

CE 
CO 


CO 
LO 


CO 

"~3 


C£ 

»— 

o 

LU 
0_ 
lO 


OO'IOI 


o 

•H 
-P 

CD 

C 

hD 

cd 

6 

CD 

CD 

fn 
.C 
P 

£ 
P 

•H 

s 
p* 
u 
■p 

o 

0) 
CO 

fn 

CD 

cd 
£} 

CO 
CO 

:o 
S     • 

CO 

C  CD 

O  P 

Sw  -H 

M  CO 


VD 


CD 

bO 
•H 
Ph 


00*001 


00*66  00*86  00-^6  00*96 

JU.ISN31NI    3AIlbl3il 


00-S6 


00'*6 


25 


o 
o 


CO 
C\J 

r- 

Ui 
LO 
O) 

II 

UJ 


CO 

cr 

CO 


CM 

r- 
o 
o 

CD 

o 

r* 

ID 

or 

i— 

CJ 

LU 
Q_ 

in 


:  O 


•oo 

I 


o 
o 


00-301     00*00T      00*86     00*96      00*fr6     00*26 

AJ.ISN3J.NI  3AIibl3y 


00*06 


+? 


T5 

C 

a 
o 

•H 
-P 
CD 
C 

hO 


O 
-P 

■p 

•H 

e 
•a 

in 

•P 

o 

0) 

ta 
Fh 

(D 
Cd 


(D 

G 

•H 


<u 


CO  r-{ 
CO    hD 

:0   C 

S  -H 

CO 

O  CD 
h  C 
H    O 


CD 
Fh 

50 
•H 


00*89 


26 


CO 

■*■ 
cm 

CO 


II 

(Jj 


CO 

ex 

QQ 


csj 
to 


Z5 

»— 
o 

UJ 
Q_ 
CO 


08-001 


a 

o 

in 

Cm    <D 

o  c 
o 

C7 

o 

CQ 

CD  T3 

■* 

o    sit 
nt   an 

o 
o 

• 

^    CD 

-P  -H 

m 

with 
grad 

o 

<& 

o 

6  H 

CO 

—* 

3    d) 

CO 

M  -H 

\ 

-p  Cm 

E 

O 

E 

CD    O 

o 

•w 

ft«H 

o 

CQ     M 

■ 

>- 

-P 

I— 

M      O 

•-^ 

0    CD 

o 

3  rH 

• 

o 

03    CD 

CD 

o 

_/ 

X3 

<J 

OLU 

CO   +5 

•H 

• 

~> 

CQ    C 

H 

CJ 

:0  CD 

S    Pn 

CD 

CD 

H 

C    <M 

hO 

o 

O  Cm 

C 

o 

m  -H 

•H 

~ 

H  T3 

CQ 

I 


-en 
i 


CO 


CD 

M 

bO 
•H 


Of  001 


00*001  09*66  03*66  08*86 

A1ISN31NI  3AIlb13y 


0*'86 


00-86 


27 


The  program's  input  is  very  similar  to  the  NBS 
unconstrained  curve  fit  program,  PARLORS  [31  and  it  may  be 
used  in  the  same  way  by  fitting  a  singlet  "site"  to  each 
line  in  the  spectrum.  The  "objectivity"  inherent  in  the 
unconstrained  mode  may  often  be  preferable  for  simple, 
resolved  spectra.   However,  the  constrained  output  is  given 
in  terms  of  the  real  parameters  of  interest  (i.e.  shift, 
magnetic  field,  quadrupole  splitting)  and  their  uncertainties, 
which  is  usually  more  convenient  than  the  positions  and  uncer- 
tainties obtained  with  the  unconstrained  fit. 

Just  as  in  PARLORS,  program  CFIT  requires  no  initial 
estimates  of  amplitudes  or  of  the  three  parabolic  background 
parameters.   The  function  to  be  fitted  is  linear  in  these 
parameters,  and  they  are  estimated  by  a  linear  fit  before  the 
nonlinear  fitting  iterations  are  begun. 

On  the  NBS  Univac  1108,  the  program  requires  313000 
words  of  core,  with  double  precision  and  limits  of  50  para- 
meters (9  sites),  and  500  data  points.   Running  times,  with- 
out compilation,  have  ranged  from  about  30  seconds  for  spectra 
like  Figure  15  to  5  minutes  for  those  like  Figure  16 .   Changes 
in  the  thickness  broadening  formulation  are  expected  to  speed 
the  program  considerably. 

(J.  C.  Travis) 

E .   Nickel-6l  Hyperfine  Interactions 

With  the  He  cryostat  described  in  Section  2B1 ,  it 
became  possible  to  measure  internal  magnetic  fields  and 
quadrupole  splitting  in  many  nickel  complexes.   Of  particular 
interest  was  the  measurement  of  internal  magnetic  fields, 
since  Mossbauer  spectroscopy  is  at  present  the  only  method 
available . 

1 .   Internal  Magnetic  Fields  in  Nickel  Compounds 
The  internal  magnetic  field  can  be  represented  by 

H.  .  =  H   +  H,  +  H-,  .  (Q) 

int    c    1    lat  v°J 

28 


where  H   is  the  core  polarization  field,  H-,  the  orbital 
c  ^  '   1 

contribution  and  H,  ,  the  crystal  lattice  field,  the  core 

polarization  field  is  about  -275  kG,  which  is  considerably 

57 
lower  than  that  in   Fe .   The  orbital  contribution  to  the 

internal  field  is  positive,  and  is  related  to  the  g-factor. 

The  observed  internal  fields  of  a  series  of  nickel  compounds 

are  listed  in  Table  3»  and  some  of  their  spectra  are  shown  in 

1  Figures  19  to  22.   The  internal  field  is  lower  for  more 

!  covalent  compounds.   In  Nilp,  the  orbital  contribution  cancels 

the  core  polarization  field.   NaNiO„,  which  was  kindly  made 

available  by  Dr.  P.  P.  Bongers  of  Phillips  of  Eindhoven, 

Netherlands,  has  nickel  in  a  formal  oxidation  state  of  3.   It 

is  a  spin  1/2  compound,  and  the  field  is  52  kG,  compared  to 

the  100  kG  field  found  for  NiO  with  spin  1.   Further  work  is 

planned  to  correlate  the  covalency  factor  with  the  internal 

magnetic  field. 

Table  3.   Internal  Magnetic  Fields  Measured  by  Nickel 
Mossbauer  Spectroscopy  for  Selected  Nickel 
Compounds 


COMPOUNDS 

LIGANDS 

SYM 

VK 

) 

Neff 

Ho 

(kG) 

NiO 

0 

6 

h 

520 

4.6 

98 

±  3 

NiFepO^ 

0 

6 

h 

858 

96 

±  3 

NiCrpO^ 

0 

6 

Td 

78 

38 

±  7 

NiSO^ 

0 

6 

°h 

15. 

5 

5.7 

24 

±  3 

NiS 

S 

6 

°h 

200 

16 

±  21 

KNiF 

F 

6 

h 

275 

4.7 

60 

±  8 

NaNi02 

0 

6 

h 

20 

52 

±  2 

NiF2 

F 

6 

°h 

80 

3-5 

99 

NiCl2 

CI 

6 

°h 

36 

3.3 

36 

±  6 

NiBr2 

Br 

6 

°h 

65 

3.0 

75 

±  3 

Nil 

I 

6 

°h 

75 

0 

±  10 

NiCl2-6H20 

CI 
0 

2 
4 

°h 

5. 

34 

77 

±  9 

Ni(HC02)2 

0 

6 

°h 

15. 

7 

41 

±  6 
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2 .   Quadrupole  Interactions  in  Planar  Compounds 
It  has  been  shown  [4]  in  the  analysis  of  the 
electric  field  gradient  (EFG),  for  nickel  compounds  that  the 
largest  quadrupole  splitting  can  be  expected  for  planar 
nickel  compounds,  which  yield  a  Y„„-- (8/7 )<qr   >d .   The 
spectrum  of  NapNi(CNK  showed  a  narrow  resonance  line. 
Further  studies  with  the  substituted  nickel  glyoximes 
gave  a  positive  AEn  for  the  nickel  dimethylglyoxime  and 
nickel  heptoximes;  and  a  negative  AEQ  for  the  nickel 
diphenylglyoximes  which  indicates  a  large  ligand  contri- 
bution to  the  EFG .   Curve  fitting  of  the  planar  compound 
spectra  is  complicated  by  the  Karayagin-Gol ' danskii  effect, 
which  changes  the  relative  intensity  of  the  triplet.   The 
intensity  is  represented  by 

I(K,L)  =  -/   Ke  "L  cos  9  d(cos  9)  (9) 

2    2       2 
where   L  =  k  (<Z  >  -  <X  >  ),  L  >  1  for  planar  compounds 

2 
K  =  2"cgs  9   for  M  =  ±1/2 

5  2 
K  =  7  ^gs  6   for  M  =  ±3/2 

v  1+COS0    _     ..      .r-/o 

K  =  -r for  M  =  ±5/2 

2 
Substituting,  and  expanding    -Lcos  0     into  a  series 

expression   gives: 

Tf  +  -,/pN_Y**2  .  (-l)n   2n+5  Tn  (10) 

I(±l/2)-2^5    -JTi   '  (2n+l)(2n+3)   L 
n 

I(  +  C5/2Y=V*1  .  iz±±n.    3n+5  Tn  (11) 

iU^J  Z^5    n!     (2n+l)(2n+3)   L 
n 

Tf  +  tr/pN=\T^2  .  (-l)n   3n+5  Tn  (12) 

K-5/2)  2v5    -in"   (2n+l)(2n+3)   L 

n 
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for  the  triplet  intensities.   Computer  curve  fitting  is  in 

progress  to  extract  the  quadrupole  splitting  and  Karyagin 
effect  from  the  nickel  glyoxime  spectrum. 

(J.  J.  Spijkerman  and  J.  C.  Travis) 

P .   Mossbauer  Spectroscopy  of  Iron  Doped  TiO^, 
1.   Introduction 
The  objective  of  this  investigation  is  to  use  Mossbauer 

spectroscopy  as  a  nondestructive  technique  for  the  deter- 

2+    ?  + 

mination  of  the  Fe   /Fe   ratio  in  iron  doped  rutile  single 

crystals  such  as  those  used  in  masers .   The  0.05  mole  %    iron 
concentration  normally  used  in  maser  crystals  was  increased 
to  about  0.5%    in  order  to  obtain  good  spectra  in  a  reasonable 
period  of  time.   An  additional  departure  from  normal  maser 
configuration  was  the  use  of  powdered  samples  to  supplement 
the  single  crystal  data. 

The  Mossbauer  spectra  of  the  powdered  specimens  yielded 
at  least  two  nonequivalent  sites  for  each  of  the  formal 
oxidation  states.   Further  experimental  studies,  including 
temperature  studies,  are  currently  in  progress  to  attempt 
to  firmly  locate  the  observed  sites  in  the  TiO„  lattice.   The 
site  identification  program  further  includes  theoretical 
lattice  sum  calculations  of  the  quadrupole  interaction, 
including  temperature  dependence  based  on  both  the  ideal 
TiOp  lattice  and  on  various  possible  configurations  including 
additional  doping  materials. 

The  complex,  multisite,  multi-oxidation  state  spectra 
obtained  in  this  study  have  been  a  primary  stimulus  in  the 
development  of  the  constrained  curve  fitting  program  dis- 
cussed elsewhere  in  this  report.   The  spectral  complexity 
also  led  to  the  use  of  powdered  samples  to  simplify  the 
intensity  relationships  in  the  spectra.   Once  the  line 
position  identifications  had  been  firmly  established  with 
the  powder  spectra,   the  single  crystal  spectra  were 

employed  to  yield  orientation  information  about  the 

35 


electric  field  gradient  tensor,  which  added  further 
evidence  for  the  site  identification. 

Previous  work  on  the  interstitial  diffusion  of  metals 
in  rutile  has  been  done  by  Johnson,  [5],  who  used  lithium 
as  a  doping  material.   A  strong  anisotropy  in  Li  diffusion 
has  been  observed  by  these  workers  and  substantiated  by 
lattice  potential  calculations  [6].   Iron  diffusion  has 
also  been  investigated  [7],  yielding  a  diffusion  constant 
estimate  for  diffusion  parallel  to  the  crystallographic 
c-axis . 

The  previous  application  of  Mossbauer  spectroscopy  to 
iron  containing  rutile  [8]  concerned  spin  relaxation  of 
internal  magnetic  fields,  and  were  evidenced  slightly  in 
some  of  the  work  reported  here. 

2 .   Experimental 

Both  single  crystal  slices  and  powders  were  doped, 
by  several  different  methods,  with  90+%  iron-57  enriched 
iron.   The  process  yielding  the  most  uniform  doping  as 
determined  by  microprobe  analysis,  was  the  embedding  of 
the  crystal  in  a  powder  mixture  of  TiO?  and  FepO~  with 
subsequent  heating.   Various  specimens  were  then  heated 
in  different  atmospheres  (Ar,  Ar+Hp,  and  Op)  for  varying 
lengths  of  time  at  temperatures  between  900  and  1000  C. 
The  latter  treatment  was  repeated  between  Mossbauer  experi- 
ments for  several  of  the  samples,  generally  with  noticeable 
effect  on  the  spectrum. 

The  Mossbauer  experiments  employed  the  previously 
described  NBS  spectrometer,  a  Janis  liquid  helium  Dewar, 

and  a  modified  version  of  the  NBS  "chicken  feeder"  liquid 
nitrogen  Dewar  [91.   The  modification  consisted  of  a  sand- 
wich of  a  heating  element  (carbon  cloth)  between  layers  of 
insulator  which  interrupted  the  copper  heat  path  between 
the  sample  mount  and  the  liquid  nitrogen  reservoir.  (See 
Section  4-B-2  this  report).   A  temperature  range  of  110  K 

36 


to  room  temperature  was  available  with  this  system.  The 
power  to  the  heating  element  was  controlled  by  a  variable 
transformer  and  the  temperature  of  the  sample  holder 
was  monitored  with  a  copper-constantan  thermocouple. 
Temperature  stability  during  an  experiment  was  checked 
either  by  frequent  thermocouple  readings  or  by  monitoring 
the  thermocouple  with  a  strip  chart  recorder. 

Normal  transmission  geometry  was  used  for  most  of  the 
experiments  although  backscattering  was  used  on  single 
crystals  to  investigate  surface  effects.   It  was  also  hoped 
that  lower  concentrations  of  iron-57  could  be  examined  with 
backscattering  geometry,  but  the  expected  improvement  in 
signal-to-background  (percent  effect)  for  a  given  concen- 
tration was  more  than  offset  by  the  reduction  in  overall 
counting  rate. 

The  geometrical  count  rate  loss  accompanying  the  low 
temperature  and  backscattering  experiments  was  more  than 
compensated  during  later  experiments  by  the  use  of  a 
Reuter-Stokes*  sealed  krypton-COp  proportional  counter, 
and  the  improved  NBS  counting  circuitry. 

3 .   Results 

The  spectrum  obtained  from  a  melt  doped  powder  con- 
taining about  1%   Pe,  shown  in  Figure  23,  closely  resembles 
the  spectrum  obtained  by  P.  Wintersteiner  [10]  from  one  of 
M.  Alam's  specimens.   Alam's  powder  sample  was  doped  to 
0.25%  Re  concentration  during  crystal  growth  in  order  to 
assure  substitutional  doping  in  the  Ti    sites.   Both 
samples  in  Figure  23  were  doped  with  FeJ  .   The  spectrum 
similarity  implies  that  even  the  sample  with  1%    concen- 
tration (Pigure  23)  reveals  the  presence  of  only  substi- 
tutional sites.   Of  course,  the  possible  existence  of 


*Por  disclaimer  of  equipment  and  materials  see  last 
paragraph  of  Introduction. 
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Figure    23 . 


Comparison  of  1.0  mole  %   oxidized  powder 
spectrum  with  0.25  mole  %   growth-doped 
spectrum  of  M.  Alam. 


interstitial  iron  ions  with  extremely  low  recoil-free 
fractions  cannot  be  ruled  out . 

The  possible  existence  of  such  low  f-factor  inter- 
stitial sites  is  reinforced  by  the  room  temperature  and 
liquid  nitrogen  temperature  spectra  of  the  same  sample 
after  reduction  and  re-oxidation,  (shown  in  Figures  24  and 
25).   The  widely  split  site  barely  resolvable  in  the  low 
temperature  spectrum,  Figure  24,  is  only  slightly  visible, 
but  not  fitted,  in  the  room  temperature  spectrum.  Figure 
25.   Further  examination  of  the  nature  of  this  site  ob- 
viously requires  experiments  at  lower  temperatures.   These 
are  planned. 

The  reduced  form  of  the  same  sample  yielded,  at  300  K, 
the  spectrum  shown  in  Figure  26.   A  reasonably  good  fit 
to  the  spectrum  was  obtained  with  the  constrained  program 
by  assuming  the  presence  of  two  wide  doublets  and  a  sing- 
let.  The  "singlet"  (probably  a  narrow  doublet)  may  be 

3  + 

attributed  to  residual  FeJ  ,  and  the  doublets  fall  within 
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?+ 

a  parameter  range  characteristic  of  Pe   .   The  least- 
squares  fit  to  the  temperature  dependent  series  with  this 
sample  are  shown  in  Figure  27.   The  line  positions  and 
amplitudes  are  shown  graphically  in  Figure  28.   The  outer 
doublet  may  be  seen  to  have  the  type  of  monotonically 

increasing  quadrupole  splitting  with  decreasing  temperature 

2  + 

that  is  expected  with  Fe   distorted  octahedral  systems,  but 

the  inner  doublet  shows  a  random  behavior  with  temperature, 
with  little  or  no  net  change  in  splitting. 

The  success  experienced  in  the  work  with  powdered 
samples  has  not  been  equaled  to  date  in  the  work  with  single 
crystals.   This  results  primarily  from  the  fact  that  quad- 
rupole doublets  are  symmetric  for  ideal  powders  but  exhibit 
an  orientation  dependence  of  the  peak  intensity  ratio  for 
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Figure    27 . 


Temperature  dependence  of  reduced  iron  doped 
rutile  powder. 
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single  crystals.   The  additional  parameter  would  provide 
useful  data  in  simple  spectra,  but  for  the  multiple  site 
spectra  encountered  in  this  study,  the  identification  of 
the  corresponding  lines  of  doublets  becomes  overly  compli- 
cated.  It  is  hoped  that  analysis  of  the  single  crystal 
data  will  be  made  possible  by  the  results  of  the  powder 
sample  studies.   The  most  consistently  occurring  doublet 
in  the  single  crystal  data  is  shown  in  Figure  29,  spectra 
(a)  through  (f),  with  powder  data,  (g)  and  (h),  shown  for 
comparison.   These  doublets  may  be  interpreted  as  substi- 
tutional  Fe   ,  although  they  are  consistently  smaller  and 
left  shifted  from  the  correspondingly  identified  site  in  the 
powder . 

A  particularly  interesting  sequence  of  single  crystal 
results  is  shown  in  Figure  30.   The  crystal  used  for  these 
spectra  was  originally  about  .8  mm  thick,  yielding  spectrum 
(c),  a  narrow  doublet  (Fe   ,  substitutional).   Because  of 
the  relatively  great  attenuation  of  the  thick  crystal, 
backscattering  was  attempted  in  the  hope  of  improving  the 
signal-to-noise  ratio.   The  unexpected  result  was  the  com- 
pletely different  spectrum  shown  in  (b),  probably  due  to 

2  + 

Fe   .  The  crystal  was  then  lapped  to  a  thickness  of  .2  mm, 

yielding  spectrum  (a),  which  contains  both  of  the  previously 
observed  doublets.   These  results  are  taken  to  imply  a 
depth  dependence  of  the  relative  concentrations  of  the 
Fe^+  and  Fe    states.   Further  controlled  experiments  of 
this  type  are  underway  to  investigate  this  phenomenon  further. 

4 .   Theoretical 

Although  oxidation  state  identification  may  be  readily 
accomplished  by  use  of  the  isomer  shift,  site  identification 
has  been  a  more  difficult  matter.   A  computer  program  for 
the  theoretical  prediction  of  quadrupole  splittings,  written 
by  Dr.  M.  J.  Bielefeld,  (See  Section  G),  is  being  employed  for 
this  purpose.  This  program  performs  lattice  sum  calculations 
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Figure  29.   Comparison  of  iron  doped  single  crystal  and 
powdered  rutile. 
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Figure  30.   Spectra  of  single  crystal  rutile  as  a  function 
of  depth. 


of  the  ligand  electric  field  gradient  tensor  and  the  potential 

2+ 

experienced  by  the  d  electrons  for  an  "Fe    ion  in  a  point 

charge  lattice.  The  program,  in  addition,  calculated  the 
d  orbital  splittings,  and  the  electric  field  gradient 
tensor  due  to  the  above  half-full  shell  Boltzman  population 
of  the  d  orbitals  for  any  desired  temperature .   The  ligand 
and  valence  contributions  are  appropriately  summed,  the 
resultant  tensor  diagonalized,  and  the  quadrupole  splitting 
calculated.   Only  a  moderately  good  fit  to  the  experimental 
data  was  obtained  by  adjusting  the  low  temperature  limit, 
by  suitable  choice  of  nuclear  quadrupole  moment  and/or 
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Sternheimer  shielding  factor  (neither  is  an  accurately 
known  quantity)  and  by  adjusting  the  d  orbital  spacings 
by  suitable  choice  of  "effective  charge"  on  the  ligands 
(a  reflection  of  convalency  effects). 

Another  type  of  fit,  a  least  squares  Boltzman  fit  based 
on  the  assumption  of  tetragonally  distorted  octahedral 
symmetry,  was  also  obtained,  and  is  shown  as  the  solid  line 
of  Figure  31-   The  d  orbital  splitting  parameter  (A™) 
is  illustrated  in  Figure  32,  and  the  Boltzman  function  for 
the  quadrupole  splitting,  neglecting  lattice  contributions, 
is  given  by 


A(T)  =  A(0)- 


1  -  exp(-AT/kT) 
1  +  2  exp(-  AT/kT) 


(13) 


where  A(0)  is  the  low  temperature  (T=0)  limit  of  the  quadrupole 
splitting.   The  fit  shown  in  Figure  31  was  obtained  for  para- 
meter values  of  AT=522  cm"  and  A(0)=2.69  mm/s . 
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Figure  31. 


Temperature  vs  splitting  of  outer  doublet  for 
iron  doped  reduced  form  of  powdered  rutile . 
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The  theoretical  calculations  corresponding  most  closely 
to  the  Boltzman  fit  were  those  using  an  effective  charge 
parameter  of  .75  (i.e.  assuming  the  point  ions  to  bear  only 
75%  of  their  normally  assumed  "valence"  charges)  for  which 

the  two  lowest  d  orbital  splittings,  A   and  A   of  Figure 

-1 
32,  were  459  and  535  cm   .   Low  temperature  limit  may  be 

normalized  to  the  least  squares  fitted  value  of  2.69  mm/s 
by  choice  of  an  abnormally  low  value  of  Q(l  -  R)  of  0.12. 
The  results  of  this  process  are  shown  as  the  open  circles 
of  Figure  31.   Neither  of  the  fits  is  particularly  impres- 
sive but  additional  data  are  being  accumulated. 

(J.  C.  Travis  and  P.  R.  Stampfl) 
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Figure    32. 


Energy  level  diagram  of  3  d  levels  used  for 
interpretation  of  data  in  Figure  31. 
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G .   Towards  a  Hartree-Fock  Crystal  Field  Treatment  of  the 

Electrons  on  Mossbauer  Nuclei 

In  order  to  include  a  crystal  field  of  point  charges 

in  the  self-consistency  scheme  of  a  Hartree-Fock  treatment 

of  electrons  about  a  Mossbauer  nucleus,  certain  basic 

integrals  must  be  evaluated.   The  following  is  one  analytical 

method  for  this  evaluation. 

The  many-electron  problem  can  be  reduced  to  parts  which 

| are  described  in  terms  of  a  single  electron  about  the 

; Mossbauer  nucleus  being  perturbed  by  a  point  charge.   A 

function  which  describes  the  spatial  character  of  the 
I 
electron  can  be  written  as 

|n,^m>  =  -ir-elm*0    <R  (14) 

where©,,    is  a  normalized  associated  legendre  function  [11] 

and  (R      n    is  a  normalized  radial  Slater  type  orbital  [12]. 

For  example,  | 2 ,1 ,0 > denotes  a  2pQ  atomic  orbital. 

-> 
The  interaction  between  the  electron  at  a  position  r 

and  the  point  charge   q   at  a  position   R   can  be  written 
as   <n,£,m| 0 | n' , I* ,m'> 

where  0  =  ^   q  _^  (15) 

|R  -  r| 

The  functions  |n,£,n>  are  centered  in  a  coordinate  system  in 
which  the  positive  Z-axis  is  directed  towards  the  point 
charge.   Functions  in  an  arbitary  coordinate  system  can  be 
readily  transformed  into  this  system  by  a  rotation  [13]. 
Since  for  a  point  charge  the  value   R-r   has  no  <p   dependence 
in  this  coordinate  system,  the  following  simplification  holds: 


«n,*.m|0|n-.t-,m->   =  S^.-  ©£>m«njJl|o|  8,.>m.  <V,*'>(16) 

where  6    ^is  the  Kronecker  delta  function. 
m,m 
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In  order  to  separate  variables  before  integration,  these 
functions  are  expressed  in  elliptical  coordinates  [14]   with 
the  Mossbauer  nucleus  at  one  focal  point  and  the  point  charge 
at  the  other  focal  point.   The  necessary  relations  for  this 

transformation  (r,0,<|)  ->  y,v,<jO  are  the  following: 

cf>    =    (j>  ;  (    0<  <j>  <  2tt  )  ; 

|r|    =   R   (y+v)/2  ;(    1<  y  <°°    )  ; 

|8-r|    =  R   (y-v)/2  ;(-1<v<1    );(17) 

cos    9    =    (l+yv)/(y+v)  ; 

sin  9  =  /(yz-l)(l-vz)/(y+v); 


R 


3 


sin  8  r2  dr  d6  =  -S-  (y2-v2)  dy  dv 


A  typical  integral,  e.g.,  ls-2s,  Is  here  worked  in  detail. 


< 


1,0,0|S|2,0,0>  -  |^(i-)*(i:)ij/jGT-^pxp[-(i-+  £)r]  (18) 


Where  £  and  £  are  the  orbital  exponents  of  the  bra  and  ket 
functions,  respectively.   After  a  change  to  eliptical 
coordinates  and  with  the  energy  and  distances  in  atomic 
units  this  expression  becomes 

<l,0,0|o|2,0,0>=r^^(Rc)^(Rr)^/dy/dv(y+vfexp[-|(R^+|^)(y+v)] 

(19) 

<1,0,0|0|2,0,0>  =  3J7S  §  (RC)*(RC')2  [A2B0+  2A1B1+AQB2]   (20) 


where  A  and  B   are  easily  evaluated  [15]  and  are  defined 
mm  J 

below. 
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For  ease  of  notation  the  following  symbols  are  defined 
for  the  integral  <n, I ,m| 0 | n" 3l' ,m'>    : 

Aj  -f  aV  J  exp  [- 1  (|i  +  ^-)v]-,  (21) 

-/'dvv^exp  [-1  (St+|tl)v]; 


B. 


[2n+l]/2 
C   =  (RC) 


[2n'+l]/2 

D   =  (RC') 


Accordingly,  the  twenty-eight  integrals  of  interest  for  a 
metal  ion  in  the  first  transition  series  are  listed  as 
Appendix  I  in  this  report. 

A  word  of  caution  should  be  mentioned  concerning  these 
integrals.   Their  numerical  evaluation  for  specific  values 
of  R,  q,  C  and  C *    can  be  readily  performed  on  a  digital 
computer  provided  the  arithmetic  is  done  in  double  precision. 

(Michael  J.  Bielefeld) 
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3.   NUCLEAR  CHEMISTRY 

A.  Introduction 
The  program  has  continued  in  its  dual  nature  in  which 

nuclear  reaction  studies  have  been  supported  by  statistical 
investigations  in  order  to  improve  the  accuracy  of  the 
resulting  conclusions  by  finding  superior  methods  for 
treating  the  experimental  data.   Generally,  both  aspects 
of  the  program  have  contributed  to  the  solution  of  problems 
arising  in  other  research  groups  at  NBS.   Examples  are  the 
gain  shift  problem  described  in  Section  C  and  the  estimation 
of  variance  in  Section  P. 

B .  The  Discovery  of  Errors  in  the  Detection  of  Trace 
Components  in  Gamma  Spectral  Analysis 

The  following  is  a  condensation  of  a  talk  given  at 
the  1968  International  Conference  on  Modern  Trends  in 
Activation  Analysis,  held  at  the  National  Bureau  of 
Standards  in  October,  1968. 

I.   Introduction 
Least  squares  solutions  to  over-determined  systems  in 
physical  science  generally  provide  parameter  estimates,  their 
standard  errors,  and  various  measures  of  the  quality  of  the 
fit  of  the  set  of  observations  to  the  model.   Such  an  approach 
has  been  applied  for  a  number  of  years  to  the  estimation  of 
individual  radionuclide  contributions  to  an  observed  composite 
gamma-ray  spectrum  [16].   Perhaps  the  most  popular  indicators 
for  misfit  of  the  observations  to  the  assumed  model  are  the 
resulting  series  of  normalized  residuals  and  the  statistic, 
X2,  which  is  derived  from  this  residual  set.   Too  large  a 
value  for  x2  is  indicative  of  trouble,  while  a  plot  of  the 
residuals  may  often  pinpoint  the  source  of  trouble.   Errors 
which  are  frequently  uncovered  in  the  process  include  miss- 
ing components,  shifts  in  peak  position  (gain,  zero),  and 
gross,  isolated  mistakes  (dropped  channels). 
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Figure  33.   Plot  of  normalized  residuals  from  the  least 

squares  analysis  of  22Na,  54Mn,  and  background 
Channel  zero  is  the  timing  channel,  and  the 
enlarged  points  refer  to  channels  containing 
systematic  errors. 


That  errors  may  not  always  be  directly  located  by  the 
residuals,  however,  is  illustrated  in  Figure  33.   This  figure 

resulted  from  the  least  squares  analysis  of  a  mixture  of 

22         54  /22 

background,    Na,  and   Mn.   The  negative  (   Na)  Deak  and 

5I1 

positive  (   Mn)  peak  in  the  residual  plot  clearly  indicate 

the  presence  of  trouble,  but  they  do  not  directly  point  to 

its  origin.   In  fact,  the  (unintentional)  mistakes  occurred 

54 
in  the  standard  spectrum  for   Mn  in  channels  62,  65,  and  67. 
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A  more  subtle  kind  of  error-detection  problem  may  arise 
when  one  or  more  components  are  present  at  "trace"  levels  - 
i.e.,  if  they  are  scarcely  detectable.   In  order  to  examine 
this  problem  we  must  first  discuss  the  question  of  detection 
limits  as  related  to  the  least  squares  approach. 

II.   Detection  Limits  in  Overdetermined  Systems 
The  general  problem  of  limits  for  decision,  detection, 
and  determination  has  been  discussed  previously  for  the 
one  component  situation  [17].   For  the  particular  case  of 
radioactivity  (Poisson  distribution)  the  "critical  level" 
was  seen  to  be, 

L   =  k  c  (22) 

C  a  o 

and  the  "detection  limit", 

hD   "  LC   +  V*  (23) 

where:  a  .a   represent  the  standard  deviation  when  the 

o  D 
component  of  interest  is  absent,  and  present  at  its  detec- 
tion limit,  respectively,  and  ka  ,k  represent  the  values  of 
the  standardized  normal  variate  corresponding  to  the 
selected  limits  for  the  errors  of  the  first  and  second 
kinds,  respectively.   In  the  limit  of  large  numbers  of 
counts  (>100),  a  ^a   ,  and  1^21^.   Thus,  for  example,  if 
the  two  types  of  error  are  each  set  at  5% >    ka=k3  =  1  •  ^5 , 
and  the  relative  standard  deviation  at  the  critical  level 
equals  .608,  while  that  at  the  detection  limit  equals  .304. 
(Note  that  Lc    is  to  be  used  for  making  a   posteriori    decisions 
regarding  detection,  whereas  Lp  is  used  to  assess  a   priori 
detectability ) . 

Equations  (22)  and  (23)  have  been  applied  to  the  measure- 
ment process  in  which  estimation  is  carried  out  by  means  of 
least  squares.   The  method  of  solution  is  iterative  in  nature 

54 


in  that  a   must  first  be  calculated  from  a  model  containing 
o 

the  component  of  interest  at  level  zero;  then  the  level  is 
adjusted  until  equation  (23)  is  satisfied.   The  results  of 
such  calculations  are  presented  in  Table  4.   The  reference 
levels  for  the  standards  and  background  used  in  deriving  the 
detection  limits  in  Table  4  were  as  follows: 

Background  28000  total  counts 

22Na  69000  total  counts 

6i|Cu  4400  total  counts 


54 


Mn  44000  total  counts 


The  detection  limit  is  thus  seen  to  depend  both  upon  the 
levels  of  interfering  components  (and  background)  and  the 

magnitude  of  the  correlation  coefficient.   If  the  model 

54       64 
includes    Mn  and   Cu,  even  at  level  zero,  the  detection 

22 
limit  for    Na  exceeds  its  value  in  the  one  component 

model.   Variation  with  level  of  interference  is  rather 
complex;  when  one  interfering  component  overwhelms  all 
others,  the  variation  approaches  the  "square  root"  rela- 
tionship discussed  in  reference  [17]. 

22 
ts  for    Na  in  various  radionuclide 


Correlation 
,   coefficient 
Component-J   Level-J      (Na,J) 


Table  4. 

Detection  li 

mixtures .a 

Number 

of 
component 

s    Lp(22Na)b 

1 

.00012° 
.0053 

2 

.0147 
.0062 
.0104 
.0260 

64 

^:cu  0  -.92 

^Mn  0  -.50 

-^Mn  1.  -.59 
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Mn  10.         -.61 


•0291        6^Cu        Ln(.284)      -.93 
54Mn         u    1.         -.62 


Number  of  observations  equals  6. 

c  Levels  refer  to  the  respective  standards  as  unity. 
Background  =  0. 
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III.   Model  Errors 

When  one  or  more  of  the  components  in  a  mixture  lies 
below  its  detection  limit  there  exists  a  significant  chance 
of  missing  such  a  component  (by  definition) .   It  therefore 
becomes  extremely  important  to  consider  the  effects  of  such 
a  nondetected  component  upon  the  estimates  of  the  remaining 
components.   If  the  model  is  correct   -  i.e.,    includes  the 
trace  component,  "A"  -  the  remaining  estimates  are  unbiased 
and  no  difficulties  arise.   On  the  other  hand,  if  the  model 
does  not  include  "A",  either  because  of  ignorance  or  because 
of  model  simplification,  the  results  are  necessarily  biased. 
The  signs  and  magnitudes  of  the  normalized  biases  are  depen- 
dent upon  the  actual  magnitude  of  "A"  and  the  respective 
correlation  coefficients.   More  seriously,  it  can  be  shown 
that,  neither  x2  nor  the  magnitudes  of  the  residuals  may 
indicate  the  presence  of  trouble  if  "A"  lies  below  its 
detection  limit. 

Errors  of  the  type  just  discussed  may  have  especial 
significance  in  the  evaluation  of  several  trace  components 
in  the  same  sample.   For  example,  in  the  Standard  Reference 
Materials  program  at  the  National  Bureau  of  Standards ,  we 
are  concerned  with  certifying  specific  materials  (including 
high-purity  materials)  with  respect  to  their  trace  constit- 
uents. It  is  quite  likely  if  an  incomplete  model  is  used 
that,  besides  the  systematic  errors  induced  in  macro- 
constituents,  trace  components  may  be  seen  when  they  are 
absent  and/or  missed  when  they  are  present.   The  probability 
of  such  mistakes  can  become  quite  large  when  large  (positive 
or  negative)  correlation  coefficients  are  involved,  as  the 
following  example  demonstrates. 

Table  5  is  based  upon  least  squares  analyses  of  a  gamma- 

64     22        S4 
ray  spectrum  constructed  from   Cu,    Na,  and   Mn  at  levels 

Lr    (.1428),  0,  and  1,  respectively.   (Levels  are  given  in 

fill 
terms  of  the  "standards"  as  stated  earlier) .   Since    Cu 

56 


/T  1 1         O  O 

Table  5.   Effects  of  model  errors  on  the  mixture    Cu,    Na , 
and  5%n. 

22Na  Errors  (p=-.93)        Mn  Errors  (p=+.46) 

Model     Systematic   Random  (a)  Lq   Systematic  Random  (a) 

Correct 

(Includes 

64Cu)       0        .0088     .0145     0       .0082 

Incorrect 
(Excludes 

64cu)     +.0137      .0032     .0052   -.0063     .0072 


was  present  at  its  critical  level  it  will  be  missed  with 

a  probability  of  50%.      It  may  be  seen  that  the  effect  of 

54 
the  incomplete  model  on  the  macro-constituent,    Mn,  is 

a  slight  reduction  in  the  standard  error  estimate  and  the 

introduction  of  a  serious,  negative,  systematic  error. 

22 
The  effect  on  the  trace  constituent    Na  -  which  is,  in 

fact  absent  -  is  profound!   A  rather  large  decrease  in  the 

standard  error  estimate  (and  critical  level)  is  accompanied 

by  a  large  positive  systematic  error  which  causes  the 

absent  component  to  be  "detected".   It  was  observed  that 

neither  the  magnitude  of  the  residuals  (extreme=-l . 4 )  nor 

the  value  of  x2  (P  =  '72)  were  indicative  of  misfit.   A 

miniature  Monte  Carlo  experiment  was  performed  with  the 

same  system  and  it  was  observed:  that  x2  exceeded  its  .975 

22 
limit  more  often  than  one  time  out  of  six;  that    Na  was 

detected  every  time;  and  that  the  most  serious  error  in 

22 

the    Na  estimate  was  equal  to  6.05  times  its  standard 

error  -  an  event  which  would  be  expected,  by  chance,  with 

_7 

a  probability  of  2x10 

The  conclusions  are,  perhaps,  obvious.   When  dealing  with 
trace  analysis  one  must  be  extremely  cautious  in  using  too 
simple  a  model,  as  the  systematic  errors  can  be  enormous, 
while  the  popular  fit-statistics  may  appear  perfectly 
satisfactory.   The  probability  of  such  pitfalls  can  be 
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assessed  rigorously,  but  this  subject  will  be  discussed 
elsewhere.   A  principal  caution  to  be  taken,  however,  is 
to  examine  the  correlation  coefficients  of  any  trace  com- 
ponents which  are  possible  but  which  are  being  considered 
for  exclusion.   If  these  coefficients  exceed  ^±0.4  it  would 
be  wise  either  to  include  such  components  in  the  model  or 
to  "orthogonalize "  by  means  of  chemistry,  coincidence 
methods,  mass  separation,  etc. 

*  *  #  * 

The  study  of  the  effects  of  missing  components  has  been 
extended  to  the  case  where  there  are  many  degrees  of  freedom 
(many  channels)  in  order  to  learn  whether,  in  this  case  also, 
components  which  are  "detectable"   but  excluded  from  the 


model  are  apt  to  yield  discernable  patterns  in  the  residuals. 

Re 
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22     137 
Results  from  a  mixture  of  three  radionuclides    Na ,     Cs, 


Mn  -  depicted  in  Figures  34  and  35  >  display  the  normalized 
residuals  resulting  from  analysis  to  two  sodium  iodide  spectra 
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Figure  34.   Normalized  residuals  from  the  analysis  of   Na , 
137cs,  5%n;  correct  model  (horizontal  lines 
are  drawn  at  ±2  a) . 
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24 
Figure  35.   Normalized  residuals  from  the  analysis  of    Na, 

137Cs,  5^Mn;  wrong  model  (horizontal  lines  are 

drawn  at  ±2o .  ) 


of  a  mixture  of  the  three  radionuclides.   The  energy  scale  in 
each  case  runs  from  zero  to  approximately  0.94  MeV.   Figure 

34  is  the  result  of  using  the  correct  model,  whereas  Figure  35 

54 
resulted  from  the  omission  of   Mn  from  the  model.   The  wrong- 
model  pattern  is  by  no  means  marked  even  though  the  actual 

54 
level  of  the    Mn  exceeded  its  detection  limit  by  50%.   Quan- 
titative evaluation  of  the  normalized  residuals  lead  to  the 
conclusion  that  the  chance  of  being  alerted  to  this  error 
in  the  model  is  only  about  0.50.   As  was  found  in  the  earlier 

study  involving  only  a  few  channels,  the  bias  in  the  remaining 

22 
components  may  be  quite  serious.   In  this  case,  the    Na 

result  contained  a  significant  positive  bias,  while  that  for 

137 

Cs  was  slightly  negative. 

(L.  A.  Currie) 
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C.   The  Significance  of  Gain  Shift 

A  temporary  discrepancy  between  the  results  of  two 
different  analytical  methods  for  the  determination  of  ^0.1% 
Co  in  Al  neutron  density  wires,  led  to  the  application  of 
least  squares  spectrum  fitting  of  neutron  activated  wires. 
Although  the  discrepancy  was  later  discovered  to  be  due 
to  an  error  in  the  Co  standard,  the  ensuing  investigation 
revealed  some  significant  aspects  of  alternative  methods 
for  fitting  gamma-ray  spectra.   Least  squares  fitting  of 
sample  #I-3B  to  Co-standard  #2,  for  example,  resulted  in 
the  residuals  shown  in  Figure  36.   (Nal  spectral  obser- 
vations were  supplied  by  D.  Becker  and  T.  Gills  of  the 
Activation  Analysis  Section. )   The  residual  pattern  is 
characteristic  of  gain  shift,  which  amounted  to  about  .kh% 
for  this  particular  case.   Least  squares  fitting  of  the 
two  cobalt  activation  standards  to  one  another  led  to  a 
ratio  of  1.036  ±.002  (±  standard  error).   Here  the  gain 
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Figure    36. 


Neutron  density  wire  spectrum  fitting  residuals 
for  sample  I-3B  vs  standard  #2  (Co.). 
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shift  was  somewhat  smaller,  amounting  to  about  0.2$.  The 
significance  of  the  shift  in  gain,  however,  may  be  seen, 
not  only  in  its  distortion  of  the  residuals,  but  also  by 
the  fact  that  the  mass  ratio  for  the  standards  was  known  to 
be  1.022Q.   Alternative  use  of  the  peak-height  method  of 
analysis  [18]  which  is  relatively  insensitive  to  gain  shifts, 
led  to  a  result  of  1.019  ±  .009.   One  especially  significant 
point  which  results  from  this  investigation  is  that  small 
gain  shifts  as  well  as  trace  components  may  remain  undetected 
while  leading  to  bias  in  the  results  of  least  squares  analyses. 

Short  and  long  term  stabilities  of  the  gamma-ray  detec- 
tion system  used  in  the  previous  studies  were  determined  as 
follows:   System  gain  was  set  at  approximately  lOkeV/channel . 
A    Cs  source  was  positioned  over  the  Nal  crystal  to  give 
a  count  rate  of   ^2x10   cpm.   A  series  of  10  minute  (live  time) 
counts  were  then  taken,  spanning  a  period  of  14  days.   Peak 
position  height,  resolution,  and  their  respective  uncertain- 
ties; peak  area,  and  the  fit  (x2/v)>  where  v=degrees  of 
freedom  of  the  points  to  an  assumed  normal-shaped  curve, 
were  all  determined  by  means  of  a  computer  program, 
"LSNORM",  described  elsewhere  (See  Section  3,  ref  9).   Long 
and  short  term  gain  drifts  as  measured  by  the  peak  channel 
position  are  presented  in  Figures  37  and  38  respectively. 
The  measurements  comprising  Figure  38  were  taken  during 
the  interval  of  time  denoted  by  the  dashed  box  of  Figure  37. 
The  estimated  long  term  drift  rate  is  0.1%/day.   Part  of 
the  drift  observed  in  Figure  38  may  be  attributed  to  a  0.5°C 
increase  in  temperature  which  occurred  during  the  10  hour 
measurement  span.   The  overall  temperature  range  for  the  14 
day  period  was  only  0.6°C. 

The  long  term  drift  rate  is  small  and  readily  compensated 
for  by  a  daily  adjustment  of  the  phototube  high  voltage. 
The  possibility  of  gain  stabilization  by  electronic  means 
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Figure  37.   Peak  channel  drift  as  a  function  of  time:  lon^ 
term  stability. 


may  be  explored  at  some  future  time.   To  see  how  precisely 
the  high  voltage  potentiometer  might  be  adjusted  a  series  of 
measurements  was  made  in  which  the  high  voltage  was  varied 
and  peak  position  for  a    Cs  source  determined.   The  results 
of  this  experiment,  presented  in  Figure  39  >  show  that  the 
adjustment  precision  is  approximately  .03%.   High  voltage 

pot  hysteresis  was  also  found  to  be  negligible. 

137 
A    Cs  source,  that  gave  an  integrated  counting  rate 

r 

of  ^2x10  ,  cpm  was  found  to  decrease  the  gain  by  about  0.h% 
relative  to  the  same  source  measured  with  poorer  geometry 

(count  rate  of  only  2.5x10   cpm.   At  more  typical  count  rates, 

4 
say  2x10   cpm,  this  effect  is  negligible. 
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Figure  38. 


Peak  channel  drift  as  a  function  of  time:  short 
term  stability. 


Possible  shapes  for  the  background  used  to  reduce  input 
gross  counts  in  "LSNORM"  are:  (1)  BG=0;  (2)  BG=constant ; 
(3)  BG  fit  to  a  straight  line  connecting  the  (equally  spaced) 
minima  on  each  side  of  the  peak  channel.   In  Figure  40  is 
plotted  peak  area  as  a  function  of  the  number  of  channels 
used  in  the  analysis.   Results  based  on  background  options 
#1  and  #3  are  given  along  with  the  computed  (x2/v)  fit. 
The  results  are  disturbing,  not  only  in  their  internal 
disparity  but  in  the  apparently  acceptable  fit  for  several 
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Figure  39.   System  gain  as  a  function  of  phototube  high 
voltage . 

of  the  data.   It  is  thus  clear  that  this  relatively  popular 
approach  for  determining  spectral  parameters  can  be  quite 
misleading.   The  reason,  of  course,  is  that  the  mathematical 
model  is  only  an  approximation,  whose  answers  are  more 
seriously  biased  for  some  parameters  (peak  area)  than  for 
others  (centroid  energy).   Our  concern  was  the  latter,  but 
the  details  have  been  given  to  stress  the  difference  between 
model  consistency  and  "truth"  in  this  particular  area. 
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Figure  40.   Peak  area  as  a  function  of  number  of  channels 
analyzed  by  "LSNORM" .   The  closed  circles 
represent  results  based  on  the  assumption: 
BG=0;  open  circles  have  had  background  sub- 
traction from  gross  data  as  described  in  test. 
Numbers  in  parenthesis  are  \2 /x>   values  . 

(L.  A.  Currie  and  H.  L.  Steinberg) 
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D .   Quantitative  Preparation  and  Analysis  of  Gas  Samples 

The  need  for  dealing  with  small,  known  amounts  of  rare 
gases  -  both  pure  and  in  mixtures  -  has  arisen  in  a  number 
of  projects.   Mixtures  of  argon  in  air  and  carbon  dioxide 
have  been  analyzed  by  means  of  gas  chromatography  in  work 
related  to  the  atmospheres  existing  in  a  glove  box  and  the 
NBS  reactor,  respectively.   Pure  samples  and  mixtures  of 
helium  and  argon  have  been  prepared  for  calibration  of  a 
static  mass  spectrometer  used  for  nuclear  reaction  studies. 
In  this  latter  case  multiple  expansions  among  calibrated 

volumes  permitted  the  preparation  of  samples  containing  10 

15 
to  10    atoms  of  rare  gas.  Discussed  below  are  some  aspects 

of  the  measurement  systems  involved  in  these  preparations. 

1 .   Argon  in  Air  and  in  Carbon  Dioxide 

The  gas  chromatograph,  employing  a  thermal  conductivity 
cell  detector  is  shown  in  Figure  4l.   This  system  will, 
hereafter,  be  referred  to  as  the  GC-TC.   The  GC  is  of  the 
dual  column  type  for  good  base-line  stability.   For  the 
experiments  discussed,  sample  constituents  were  separated 
on  a  10  ft  long  molecular  sieve  column;  oven  and  TC  Temper- 
atures were  set  at  50°C  and  60°C,  respectively;  a  TC  bridge 
current  of  230  mA  was  used;  and  the  carrier  gas  was  helium 
at  a  flow  rate  of  about  50  ml/min . 

A  glove  box,  Figure  42,  has  been  modified  to  permit 
the  establishment  within  it,  of  high  purity,  inert  atmospheres 
Li  samples  relatively  free  of  nitride  and  oxide  "skins"  have 
been  prepared  in  this  box.   The  atmosphere  within  the  box  can 
be  sampled  through  a  septum-port  located  between  the  two 
gloved  ports. 

The  response  function  of  the  GC-TC  for  0?,  N?  and  Ar 
as  a  function  of  injected  sample  volume  in  \ii    is  plotted 
in  Figure  43.   Typical  pure  Ar  and  air  spectra  are  shown  in 
Figures  44  and  45  respectively.   As  can  be  seen  0„  and 

Ar  have  similar  retention  times.   Nevertheless,  with  Figure 
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Figure  4l.   Gas  chromatograph  with  thermal  conductivity 
detector . 


43  in  hand,  the  percentage  of  air  in  Ar  can  be  deduced  from 
the  sample  volume  and  the  height  of  either  0?  +  Ar  or  Np  peaks; 
spectra  of  two  45.5  y&  samples  drawn  from  the  glove  box  at 
different  times  during  a  flushing  operation  are  presented  in 
Figures  46  and  47. 

A  (visual)  comparison  of  darkening  rates  for  freshly 
exposed  Li  surfaces  vs  composition  of  the  glove  box  atmosphere 
has  shown  that  a  proportion  of  Ar:air  of  at  least  20:1  is 
required  to  avoid  rapid  nitride  formation. 
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Figure  42.   Glove  box  setup  for  lithium  sample  preparation. 


The  carbon  dioxide  blanket  of  the  National  Bureau  of 
Standards  reactor  was  suspected  of  containing  impermissibly 
high  content  of  Ar,  (i.e.    Ar  present  in  the  blanket  as  a 
result  of  a  leak  from  the  surrounding  atmosphere  will  react 
with  thermal  neutrons  to  form  the  1.8  h  activity  -   Ar . ) 
Our  group  was  asked  to  determine  the  level  of  Ar  in  the 
COp  "blanket"  before  and  after  several  attempts  were 
made  to  eliminate  the  leak. 

The  GC-TC  was  calibrated  for  a  range  of  air-CO„  mixtures 
from  0-100%,  using  the  operating  conditions  already  given. 
The  system  was  not  observed  to  respond  to  C0„  so  that  the 
resulting  spectra  obtained  are  in  accord  with  those  of  air 
for  equivalent  volumes.   All  C0„  blanket  samples  measured 
indicated  air  percentages  of  from  15%  to  30%. 
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Figure  43.   Thermal  conductivity  cell  response  to  0?,  Np 
and  Ar  as  a  function  of  sample  volume. 

2 .   Vacuum  System  Used  to  Prepare  Ar ,  He  and  Li  Samples 
A  portion  of  the  vacuum  system  in  which  pure  and  mixed 
samples  of  Ar  and  He  were  prepared  is  shown  in  Figure  48. 
It  is  similar  to  Figure  24  in  Ref.  9  except  for  the  following 
A  lmV/full  scale,  dual  trace,  strip  chart  recorder  has  been 
incorporated  into  the  system.   Simultaneous  monitoring  of 
pressure,  (as  measured  with  a  Pirani  gauge),  and  sample 
temperature,  (as  measured  with  a  chromel-alumel  thermocouple) 
as  functions  of  time  are  now  possible.   Calibration  curves 
for  each  of  these  parameters  are  shown  in  Figures  49  and 
50.   (Note  that  1  torr=133-32  N/m2 . ) 
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Figure  48.   High  vacuum  system  with  dual-trace  recorder  for 
simultaneous  monitoring  of  sample  temperature 
and  pressure. 


For  clarity  of  discussion  a  topological  diagram  of  the 
vacuum  system  is  given  in  Figure  51.   In  it  V   represents 
that  part  of  the  vacuum  line  to  which  the  sample  is  connected, 
but  does  not  include  the  sample  volume  per  se .   The  sample 
can  be  closed  off  from  V  by  a  stopcock,  Zj,,  also  not  shown. 
Included  in  V   is  a  precise  aneroid-type  pressure  gauge  with 
a  measurement  capability  of  from  1  to  800  torr.   V   is  that 
part  of  the  vacuum  line  to  which  a  vacuum  pump  is  connected. 
The  Pirani  gauge  samples  P  . 
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Figure  h9 .      Recorder  response  as  a  function  of  vacuum  system 
pressure . 
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Figure    50. 


Recorder  response   as    a   function   of   sample   oven 
temperature . 
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Figure  51.   Topological  diagram  of  high  vacuum  system. 


Critical  factors  in  the  preparation  of  picomole 
quantities  of  gas  samples  are: 

(1)  degassing  of  samples  before  filling,  and  potential 
behavior  of  these  filled  samples  in  the  mass  spectrometer. 
A  sample  into  which  a  known  quantity  of  gas  has  been  intro- 
duced is  said  to  be  filled.   Degassing  is  accomplished  by 
heating  the  copper  capsule  (with  or  without  Li)  to  about 
300°C  under  high  vacuum.   To  test  the  efficacy  of  this  pro- 
cedure and  to,  in  effect,  determine  the  probable  behavior 

of  the  filled  samples  in  the  mass  spectrometer  several 
Cu-capsules,  containing  *\»  0 .  lg  Li  each,  were  degassed  and 
then  filled  with  He  and/or  Ar  at  pressures  of  about  100  torr. 
Each  sample  was  then  cooled  and  later  evacuated.   With  Z. 
closed  but  Z1  open  the  sample  temperature  was  gradually 
increased  to  300°C.   A  typical  run,  as  observed  by  the 
recorder,  is  shown  in  Figure  52.   Note  that  the  break- 
point in  the  pressure  curve  came  at  the  Li  melting  point 
(T=179°C). 
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Figure  52.  Pressure  and  temperature  response  curves  obtained  by  melting 
a  (O.lg)  LI  sample  In  which  Ar  at  a  pressure  of  100  torr  has  been  "frozen" 
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(2)  Determination  of  expansion  ratios,  r,  of  the  form: 
r=V  /V  +V  ) .   Neither  the  aneroid-type  nor  the  pirani 

&  lD  \^ 

pressure  gauges  are  satisfactory  at  the  pressures  required, 

_o        _o 

(i.e.  10    to  10    torr),  nevertheless,  a  method  of  obtaining 

low,  but  accurately  determined  pressures  was  devised.   Using 

the  aneroid-type  gauge  in  its  sensitive  range  the  ratio,  r, 

was  determined  for  two  V  volumes   V   and  V   ,  differing 

s  S]_      S2 

by  the  addition  of  two  bulb's  to  V  .   These  bulbs  can  be 
J  s 

closed  off  from  the  remainder  of  V  by  stopcocks.   Ideal 

gas  behavior  and  isothermal  conditions  are  assumed. 

The  calibration  of  the  volume  ratio,  r,  was  performed 

as  follows.   (Note  that  no  sample  need  be  present  for  this 

calibration.)   Volume  V  +V  +V  was  filled  with  either  He  or 

s   e   g 

Ar  gas  to  a  pressure  of  about  800  torr  (=P  )  as  measured 

by  the  aneroid-type  pressure  gauge.   Z1  was  then  closed  and 

V  and  V   evacuated.   Zn  was  then  closed  and  V  expanded 
eg  3  s   F 

into  V  +V  .   After  equilibrium  Z,  was  closed,  Z?  opened 

and  V  expanded  into  V  +V  .   The  final  pressure  P   is  then 
e  eg  ^         g 

measured  and  the  two  pressures  P   and  P  are  related  as 

s      g 


P  ,  /P  =  r   =  (V./V  .  )(V  /V    )  (?4) 

e+g   s    o    v  s   s+e    e   e+g  k^^) 

with  V  at  some  measured  pressure  P'   Z„  is  closed  and  V 
g  g    2  s 

and  V  are  evacuated.   Then  Zp  is  opened,  V  expanded  into 

V  +V  ,  and  the  final  pressure  P'    recorded.   P   and  P 
g   e5  v  g+e  g      g+e 

are  related  as 

p  '/P '    e  r   =  V   /V  roc  > 

g7  e+g    a    e+g7  g  (25) 

Equations  24  and  25  may  be  combined  to  determine  r  as  a 

function  of  r   and  r  ,  namely 
o      a 

V/V^   Er  =  r(l  +  — =kp)  (26) 

s   s+e        o     r  -1'  \<^^j 
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Quantitative  filling  of  samples  then  proceeds  as  follows. 

First  the  sample  is  degassed,  the  entire  volume  V=V  +V  +V 

&      s  s   g   e 

having  been  evacuated  in  this  process.   The  sample  is  then 
isolated  by  closing  Zu.      A  quantity  of  pure  gas  is  introduced 
into  V  and  the  gas  pressure  P  ,  determined  by  means  of  the 
aneroid-type  gauge.   Z?  is  then  closed  and  by  m  successive 
evacuations  of  V   (with  Z-,  closed  and  Z^  opened)  followed 

by  equilibrations  of  V  and  V   (with  Z-.  opened  and  Z-  closed) 

n 
P   is  reduced  to  some  final  pressure  Pf,  by  a  factor  r  . 

With  Z-,  closed  and  a  sample  capsule  present  Zu    is  opened  and 

the  capsule  filled  to  a  pressure  P'  readily  determined  from 

P„,  V  and  the  capsule  volume. 
f '   s  ^ 

Precise  mixtures  of  Ar  and  He  are  obtained  in  V 

(and  then  the  sample)  by  confining  each,  in  one  of  the  bulb 

appendages  of  V  ,  at  a  known  pressure  and  then  opening  each 

to  V  after  V  had  been  evacuated, 
s        s 

(H.  L.  Steinberg  and  L.A.  Currie) 


E.   Photonuclear  Yield  Analysis 

Copper  discs  have  been  irradiated  with  bremsstrahlung  from 
the  NBS  synchrotron.   Yields  were  obtained,  at  20  endpoint  ener- 
gies, X.,  from  19.5  MeV  to  166  MeV  for  the  following  reactions: 
^CuCyA);  63Cu  (Y,n);  63Cu(y,2n);  and  63Cu(y,3n).   This 
experiment  was  the  basis  for  one  of  the  authors  (HLS)  M.S. 

Thesis  [19]. 

Several  methods  exist  for  unfolding  yield  data  to  produce 
either  a  histograph-shaped  excitation  function  or  integrated 
cross  sections  [20].   As  discussed  in  Reference  9  existing 
methods  of  analysis  magnify  yield  uncertainties  10-  to  100- 
fold.   Alternate  methods  of  analysis  have  been  explored  in 
preparation  for  future  experimentation  by  this  group.   In 
particular  we  seek  techniques  that  will  give  definitive  or 
limiting  information  about  the  "tail"  region,  (i.e.  the  region 
E  £  30  MeV)  . 
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All  techniques  used,  both  old  and  new,  involve  a  least 
squares  (or  simultaneous  equation)  reduction  of  yield  data. 
In  one  technique  explored  the  shapes  of  the  (y,n)  and  (y,2n) 
excitation  functions  in  the  giant  resonance  (GR)  region, 
(E,,  E  30  MeV)  are  either  taken  bodily,  from  previous  results, 
or  are  assumed  Lorentzian  [21].   Above  30  MeV,  where  the 
literature  is  sparse,  the  following  shapes  have  been  used: 

(1)  c(E)  =  0; 

(2)  a(E)  =  an,    ,; 

D(y,n) ' 

E2r2 

(3)  a(E)  =  Lorentz  Shape: 

E2T2+(E2-E2)2 

(4 )  a   =  Const . ;  and 

(5)  a  combination  of  (2)  and  (3). 

If  the  GR  is  assumed  Lorentz-shaped  the  yield  equation 
can  be  written  as 

30  X. 

Y(V  =  a0/N(Xi'E)aL(E>Eo'r)dE  +  a0/1N(X.,E)T(E)dE   (27) 
Eth  3° 

=  Vi(  Wr)   +  V2(V 

By  means  of  a  linear,  weighted,  least  squares  program, 

written  for  a  time-sharing  computer,  the  amplitudes,  a   and 

a , ,  of  assumed  GR  and  tail  shapes  may  be  determined.   Our 

63cu(y,n)  yields  with  error  bars  representing  ±1  standard 

error  are  plotted  in  Figure  53.   For  comparison  a  composite 

curve  comprised  of  the  actual  a(E)  values  of  Sund,  et  al .  [22], 

to  25  MeV  and  Fultz,  et  al .  [23],  from  26  to  28  MeV,  is  given. 

If  we  assume  no  tail  (shape  number  1),  the  dashed  part  of 

the  yield  curve  results.   Although  this  curve  is  consistent 
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with  our  data  below  about  40  MeV,  our  yields  appear  to  depart 
from  this  curve  above  HO   MeV.   This  suggests  the  existence 
of  non-negligible  tail.   Our  data,  when  fit  with  this  shape, 
(i.e.  Sund  +  Fultz  +  tail  shape  number  1),  result  in  the 
normalized  residuals  depicted  in  Curve  A  of  Figure  53.   If  a 
"quasideuteron"-type  tail  (shape  number  2)  is  used,  the 
residuals  plotted  as  curve  B  result.   It  is  evident  that  by 
assuming  a  "quasideuteron"  tail  the  fit  is  improved  considerably 
This  does  not  necessarily  imply  the  correctness  of  this  shape 
so  much  as  the  need  for  additional  cross  section  above  30  MeV 
to  satisfy  our  high  energy  yield  data.   A  similar  need  for 

additional  cross  section  above  30  MeV   also  resulted  from 

65 
analysis  of  our   ^Cu(Y,n)  yield  data. 

If  E   and  V   are  not  known  a  Taylor  expansion  can  be 

performed  about  these  quantities,  thus 

Y(X.)(n)  =  c  I,(X. ,E(n),r(n))  +  c  (3I-/3E  )5E(n) 
i        olio'         o    1    o    o 

s  (n)  (28) 

+  aoOl1/9D6rln;  +  atI2(X±) 

A  nonlinear,  iterative,  weighted  least  squares  option 

exists  for  the  program  discussed.   Input  guesses  E 
/  \  o 

and  T     for  E   and  r  are  made.   Each  successive  least 

O  /       ,  V 

squares  iteration  results  in  new  estimates:  E      = 

E(n)      (n)    (n+1)  .  r(n)      (n)    (n+1)  a°d   fn+1) 

o        o  '  '   o  t 

If  either,  or  both,  test  ratios  \6    E  ^n)/  E  ^ | , | 6T (n)/r (n) | 

are  less  than  some  preset  value,  (i.e.  the  value  .001  has 

been  used  to  date),  another  iteration  is  performed  using  the 

updated  values  for  E   and  r.  In  addition  to  computing  a  , 

the  optional  parameters  E  ,  r,  and  a,   ,  and  their  respective 

uncertainties,  yields  based  on  these  "least  squares", 

parameters  are  formed  at  each  X..   Residuals,  normalized 

residuals,  average  cross  sections  for  each  interval, 

A  =  X..-.-X.  and  the  integrated  cross  section  at  each  X. 
l+l   i  to  i 
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Figure  53.   Experimental  yields  ±  1  standard  error  for  the 
reaction  °3cu(y,n)  for  comparison:  the  solid 
curve,  based  on  the  excitation  function  for  this 
reaction  due  to  Sund  [22]  to  25  MeV  and,  from 
25  to  28  MeV,  Fultz  [23].   If  cr(E)  is  assumed 
zero  above  28  MeV  the  dashed  curve  results . 
Residuals  obtained  when  our  yield  data  (to  166 
MeV)  were  fit  with  the  above  shape  are  plotted 
in  Curve  A.   If  a  "quasideuteron"  tail  is 
assumed  above  30  MeV  least-squares  fit  to  our 
data  gives  the  residuals  of  Curve  B. 
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are  determined.   Goodness  of  fit  quantities:  x2/v  and  the 
serial  correlation  coefficient  for  the  set  of  data  are  also 
part  of  the  program  output . 

(L.  A.  Currie  and  H.  L.  Steinberg) 

F .   Indirect  Variance  Estimation  in  Chemical  Systems 

A  method  for  assessing  the  variability  of  a  chemical 
component  has  been  studied  for  the  case  where  the  component 
cannot  be  directly  measured,  but  where  information  is  derived 
from  observations  on  a  second,  combining  component.   The  method 
is  based  upon  a  simple  physical  model  incorporating  stoichio- 
metry  and  the  observation  of  the  variability  of  the  second 
related  component  and/or  the  ratio  of  the  two.   A  principal 
result  of  the  investigation  is  an  expression  giving  limits 
for  the  relative  variance  of  component-A,  <|>j?: 

l*A/B  "  *Bl  1   *l<    ♦1/B  +  *B  (29) 

where  ^aad  an^  <J>R  represent  the  relative  standard  deviations 
of  the  ratio  of  components-A  and  -B  and  of  component-B, 
respectively. 

The  problem  arose  because  of  the  need  to  estimate  the 
heterogeneity  in  the  calcium  content  of  chemical  microstandards , 
but  in  which  it  was  not  possible  to  carry  out  direct  obser- 
vations of  the  calcium  variability  [24].   The  problem  was 
solved  by  examining  the  variability  of  the  ratio,  Ca/S,  and 
that  of  the  sulfur  (-S0~),  to  which  the  calcium  was  bound. 

Derivation  of  the  inequality  [24]  was  based  upon  the 
Taylor  expansion  for  error-propagation:  [25] 

*A/B  "  *A  +  *B  "  2PAB*A*B  (30) 

where  p.R   represents  the  correlation  coefficient  between 
the  variations  of  components-A  and  -B.   The  goal  of  the 
indirect  approach  is  to  estimate,  or  set  limits  for, 
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(f>.  -  which  is  presumed  immeasurable  -  based  upon  obser- 
vation of  (j)„  and/or  <{>.  ,Ronly.  (We  assume  no  information 
regarding  the  correlation  coefficient.) 

In  order  to  derive  limits  for  <J>.  using  Eq .  (30)  a 
stoichiometric  model  was  assumed,  in  which  the  two  chemical 
species  were  assumed  to  have  a  common  source  of  variation 
(fluctuations  in  the  amount  of  the  chemically  bound  substance) 
plus  independent  sources  (fluctuations  characteristic  of 
excess  amounts  of  the  components).   Based  upon  this  model, 
one  can  establish  limits  for  the  ratio,  <J>./(}>  ,  as  given  in 

.H.    ID 

Equation  ( 31 ) , 


'AB  ±  V*B  -  1/PAB  (3D 


assuming  the  mixture  of  components  to  be  stoichiometric,  on 
the  average.   Alternative  models  which  might  be  proposed 
include:  (a)  fixed  amount  of  chemically-bound  A  and  B  plus 
varying  excess  amounts  of  either  or  both,  and  (b)  a  fixed 
amount  of  B,  but  varying  fractions  available  for  chemical 
binding.   Both  of  these  latter  models  can  be  considered 
special  cases  of  the  original  model,  in  which  p.  =0.   The 
first  alternative  eliminates  the  last  term  from  Equation 
(30),  and  the  second,  the  last  two  terms. 

Equations  (30)  and  (31)  permit  the  setting  of  limits 
for  the  heterogeneity  in  component-A  as  a  function  of  p„R 
based  upon  knowledge  of  <f>R  and/or  $a/tv   ^he  limits  are 
represented  graphically  in  Figure  5^>  which  also  indicates 
some  concomitant  restrictions  on  <j>R  and/or  4>a/r- 

If  one  has  available  chemical  or  physical  information 
which  may  be  used  to  restrict  the  range  of  p  then,  of  course 
the  curves  in  Figure  5^  may  be  used  to  derive  <J>.  limits. 
which  are  correspondingly  narrower  based  upon  this  additional 
constraint . 
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Extending  the  study  to  the  case  where  cf>R  and  <J>«/R 
(but  not  PAR)  are  known  leads  to  the  inequality  given  in 
Equation  29.    One  conclusion,  which  will  be  dealt  with  in 


<f> 


A/B 


is 


more  detail  elsewhere  is  that  if  the  ratio, 

-1/2  ^R 

less  than  5     ,  <K  is  quite  well  approximated  (within 

l  /? 
about  10$)  by  <|>R,  whereas,  if  the  ratio  exceeds  5    , 


B 


* 


* 


A    rA/B 

This  study  has  been  applied  to  Ca  variation  in  chemical 

microbeads,  as  previously  stated,  in  which  the  variability 

of  sulfur  and  of  the  ratio  Ca/S  were  employed  to  deduce 

limits  for  that  due  to  Ca,  and  in  which  no  knowledge  of 

the  correlation  coefficient,  on      Q,  was  assumed.   The  lower 

limits  for  <j>p   as  reported  in  ref.  [24].  was  found  to  be 
o  a 

zero,  while  the  upper  limit  (Equation  29)   was  7%. 

(L.  A.  Currie) 
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G .   High  Sensitivity  (Radio) Chemistry 

Unique  opportunities  and  problems  associated  with  the 
detection  or  determination  of  fewer  than  one  billion  atoms 
of  radioactive  or  stable  nuclides  have  recently  presented 
themselves  as  a  result  of  (1)  a  study  for  the  ICRU  on 
The  Measurement  of  Low  Level  Radioactivity  [26],  and 
(2)  research  in  nuclear  chemistry  and  activation  analysis 
in  which  the  probability  of  nuclear  reaction,  leading  to 
a  given  residual  nucleus  is  extremely  small.   The  small 
reaction  probability  generally  results  from  a  small  nuclear 
cross  section  in  the  case  of  nuclear  chemistry,  or  small 
(trace)  amounts  or  concentrations  of  target  nuclei  in  the 
case  of  activation  analysis. 

In  the  absence  of  interference  or  contamination  a 
number  of  the  methods  used  in  this  laboratory  -  -  for 
example,  direct  measurement  of  long-lived  radionuclides, 
and  measurement  of  stable  nuclides  by  means  of  nuclear 
activation  or  static  mass  spectrometry  -  -   permit  the 

o 
detection  of  femtomole  (^6x10   molecules)  quantities.   The 

fundamental  limitations  in  this  case  include  detector 
background,  detection  efficiency  and  chemical  or  physical 
yield,  and  statistical  fluctuations  in  the  signal,  itself. 
Although  extremely  powerful  means  exist  for  detecting 
pure  materials,  reagent  and/or  sample  contamination  fre- 
quently fix  the  actual  detection  limits  at  levels  far  in 
excess  of  those  which  characterize  the  instrument  capa- 
bilities.  To  cite  an  extreme  case,  the  sea-level  flux  of 
secondary  neutrons  originating  from  cosmic  ray  interactions 
in  the  atmosphere  is  sufficient  to  produce  a  steady  state 
concentration  of  32P  in  HC1  of  about  7000  atoms/kg.  (Note 
that  the  conventional  6  lb  bottle  of  concentrated  hydro- 
chloric acid  contains  about  1  kg  of  HC1 . )   Beta  decay  of 
the  14. 3  day  J  P  in  turns  leads  to  a  linear  increase  of 

op 

(stable)  J   S  which  is  equivalent  to  an  impurity  concentration 
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of  7-^  ppT/10   years.   Although  such  an  impurity  level  of 
stable  sulfur  is  trivial,  even  for  trace  analyses,  the 
| impurity  level  of  radioactive  phosphorus  might  cause 
difficulty  in   specially  low-level  radiochemical  procedures 
in  which  quite  large  amounts  of  HC1  are  employed  as  a 
reagent  [27] . 

In  addition  to  chemical  (or  radiochemical)  Interference, 
two  forms  of  contamination  are  particular  noxious  [28], 

The  first,  isotopic  contamination,  is  illustrated  by  the 

85 
production  of    Kr  in  nuclear  fission  such  that  the  specific 

radioactivity  of  atmospheric  krypton  has  been  raised  to  about 

140  pCl/mmole.   This  contamination  has  posed  a  significant 

On 

problem  for  the  detection  of  cosmic  ray  produced   Kr , 

whose  atom  concentration  in  atmospheric  krypton  is  one 

85 
million  times  smaller  than  that  of    Kr  [29].   A  still 

more  subtle  form  of  contamination  is  isonuclidic .   Here, 
the  reagents  or  samples  are  contaminated  with  the  very 
nuclide  sought,  and  therefore,  neither  physical  nor  chemical 
methods  may  be  applied  for  separation.    H  serves  as  one  of 
the  most  characteristic  examples.   In  an  international  com- 
parison of  low  level  tritium  samples,  "Ostlund  reported 
measurements  in  which  mine  water  used  in  the  preparation 
of  the  lowest  concentration  HTO  samples  was  shown  to  contain 
a  tritium  concentration  (T/H)  of  less  than  10~    [30]. 
Due  to  the  proximity  of  nuclear  reactors,  however,  one  of 
the  participating  laboratories  found  levels  In  the  blank 
solution  some  ten  to  twenty  times  higher. 

When  one  is  faced  with  forms  of  interference,  other  than 
isonuclidic,  the  detection  limit  can  often  be  substantially 
reduced  [31]  by  chemical  or  physical  separations  according 
to  Z  (atomic  number)  or  A  (mass  number)  or  both.   Figure  55 
illustrates  a  general  scheme  for  sample  purification  in 
order  to  concentrate  the  radionuclide  of  interest  as  well 
as  to  separate  it  from  interfering  radionuclides.   The 
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Figure  55.   General  scheme  for  low  level  radiochemical 
separations  and  identification. 


four  stages  of  separation  shown  vertically  (bulk,  chemical 
group,  chemical  element  (Z),  and  mass  (A)  which  may  serve 
as  increasingly  specific  interference-reduction  alternatives 
to  non-destructive  analysis,  which  is  represented  by  the 
dashed,  horizontal  line.   Two  modes  of  "ZA"  discrimination 
are  represented  by  horizontal,  solid  lines.   Both  modes 
depend  upon  nuclear  characteristics:   the  first,  upon  the 
nuclear  decay  scheme;  and  the  second,  upon  the  unique 

parent-daughter  chemical  relationships  resulting  from 
radioactive  decay.   The  particular  radiochemical  analysis 
represented  in  the  figure,  Ce-l44  in  spinach,  has  been  more 
fully  discussed  in  the  literature  [32]. 

Mass  separation  is  highly  desirable  when  one  encounters 
isotopic  contamination.  Figure  56  illustrates  schematically 
the  approach  taken  in  this  laboratory  for  the  physical  sep- 
aration of  1-131  from  1-125  which  was  present  in  the  original 
sample  as  a  diluting  radioisotope  [33]. 
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Figure  56.   Mass  separator  experiment. 

In  order  to  give  prespective  to  the  enormous  sensitivity 
which  may  be  achieved  by  a  judicious  combination  of  the 
separation/identification  methods  depicted  in  Figure  55, 
summary  of  a  characteristic  series  of  experiments  involving 
the  detection  of  extremely  small  amounts  of  material  (10   to 
10    atoms)  has  been  collected  in  table  6.   The  experiments 
have  been  listed  in  order  of  decreasing  element  concentration, 
and  the  principal  means  (in  addition  to  bulk  concentration) 
for  purification  have  been  included.   The  first  two  experi- 
ments refer  to  nuclear  reaction  studies  underway  at  NBS;  the 
effects  of  cosmic  rays  on  HC1  have  already  been  discussed; 
and  the  last  three  experiments  represent  large  scale  studies 
related  to  cosmic  ray  interactions  in  the  vicinity  of  the 
earth.  Note  that  the  concentration  of    Ar  which  may  be 
detected  in  CCIl  as  shown  in  last  row  of  Table  6  is  less 
than  one  atom  per  80  moles  (CI)!   It  should  be  clear  from 
the  table  that  in  the  most  highly  sensitive  experiments 
even  so  rare  a  process  as  cosmic  ray  contamination  of  the 
sample  or  the  reagents  indeed  may  not  be  negligible. 

(L.  A.  Currie) 
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h.       NUCLEAR  INSTRUMENTATION 

A .  Introduction 

This  activity  has  been  transferred  from  this  section  as 
of  July  1,  1969.   A  new  group  called  the  "Technical  Support 
Group"  will  be  formed  with  staff  from  this  section  and  section 
no.  8.   To  some  degree  this  arrangement  has  been  operating 
already.   For  example,  Mr.  Shideler  has  been  detailed  from 
this  section  to  section  6  for  the  past  two  years. 

The  new  Technical  Support  Group  will  be  staff  support 
to  the  Chief  of  the  Analytical  Chemistry  Division.   Its 
members  are:   B.  Bowman,  F.  A.  Lundgren,  J.  Matwey,  T.M. 
McCleod,  F.C.  Ruegg  and  R.  W.  Shideler.   Further  information 
on  the  function  of  this  group  will  appear  in  a  new  progress 
report  next  year. 

B.  Optical  Drive 

1 .   Introduction 

A  Mossbauer  Spectrometer  of  the  constant  acceleration 
type  which  uses  an  interferometer  to  sense  and  control  the 
velocity  with  high  accuracy  and  precision  has  been  built 
at  the  National  Bureau  of  Standards.   The  spectrometer 
consists  of  five  component  parts:   (1)  the  drive  electronics, 
(2)  the  optical  interferometer,  (3)  the  frequency  ratio 
detector,  (4)  frequency  synthesizer  and  its  control  logic, 
and  (5)  the  drive  transducer,  [9]. 

In  figure  57  a  block  diagram  of  the  system  is  shown. 
The  spectrometer's  velocity  input  signal  is  derived  from 
the  multichannel  scaler's  address  analog  [39]  which  is 
a  sawtooth  voltage  waveform  vs  time.   The  drive  electronics 
amplifies  this  signal  and  compares  it  with  two  feedback 
signals,  one  from  the  linear  velocity  transducer  and  one 
from  the  optical  system.   The  feedback  from  the  optical 
system  is  derived  from  the  fringe  frequency  compared  to 
the  frequency  of  the  digitally  controlled  frequency  synthe- 
sizer.  When  the  two  frequencies  are  unequal  a  correction 
signal  is  applied  to  the  integrator  which  makes  a  continually 
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Figure  57.   Block  diagram  of  the  optical  drive. 

larger  correction  until  the  correct  velocity  is  achieved. 
The  fringes  or  the  gammas  coming  from  the  source  can  now  be 
counted.   The  fringe  count  will  be  proportional  to  the  velocity 
in  each  channel  and  the  velocity  can  then  be  calculated  from 
the  fundamental  values  of  the  wave  length  of  the  inter- 
ferometer and  time  with  a  high  degree  of  accuracy. 
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2 .   Drive  Electronics 

The  drive  electronics  are  an  improved  version  of  the 
electronics  as  published  before  [9,  ^0,  4l],  and  the  new 
circuits  will  be  described  in  detail  here.   The  electronics 
can  be  divided  into  three  parts:  (1)  drive  amplifier, 
(2)  photocell  trigger  circuit  and  (3)  logic  circuit.   The 
drive  electronics  are  designed  to  control  our  new  transducer 
[9].   The  electronics  package  is  designed  to  be  assembled 
into  a  triple  width  Nuclear  Instrument  Module,  and  it  derives 
all  of  its  operating  voltages  from  a  NIM  standard  bin  and 
power  supply.  (See  figure  58). 

The  drive  electronics  will  operate  in  a  nonoptical 
Mossbauer  spectrometer  system  by  eliminating  all  of  the 
drive  components  which  are  related  to  the  optics  and  the 
optical  feedback  path. 

a.   Drive  Amplifier 

The  schematic  of  the  drive  amplifier  is  shown  in 
Figure  59.   The  analog  of  the  address  scaler  is  brought 
in  on  Pin  21  of  the  drive  amplifier  circuit   board.   It 
is  summed  with  a  DC  (direct  current)  level  from  potentio- 
meter, R20,  so  that  it  will  be  symmetrical  with  respect 
to  ground.   The  ramp  is  then  summed  with  the  correction 
signal  produced  by  the  optical  interferometer,  and  its 
amplitude  is  adjusted  by  R22.   The  ramp  signal  is  now 
compared  with  the  output  of  the  linear  velocity  transducer 
by  operational  amplifier  number  one.   Amplifier  nos .  2  and 
3  amplify  the  signal  and  control  the  frequency  response 
so  that  it  can  be  applied  to  the  symmetrical  power  amplifier. 
The  power  amplifier  consists  of  four  transistors:  MD985  a 
dual  package  and  STC  1080  and  2N3171  power  transistors 
mounted  on  the  chasis  for  the  purpose  of  heat  dissipation. 
The  power  amplifier  has  a  voltage  gain  of  1. 
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Figure    58.      Photograph   of   Doppler  drive    electronics. 
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b .   Photocell  Trigger  Circuit 

The  purpose  of  the  photocell  trigger  circuit  is  to 
shape  the  signal  arriving  from  the  photocell  position  sensor 
mounted  on  the  transducer,  to  generate  a  fast  rise-time 
pulse  to  control  some  of  the  logic  circuits,  and  to  gene- 
rate an  address  scaler  reset  within  the  multichannel  scaler. 
The  circuit  also  has  a  bistable  multivibrator  to  generate 
the  input  signal  when  the  spectrometer  is  operating  in 
the  constant  velocity  mode. 

A  schematic  of  the  photocell  trigger  circuit  is 
shown  in  Figure  60.   The  photocell  signal  enters  the  circuit 
board  on  pin  5.   When  the  slit  is  between  the  lamp  and  the 
photocell  the  light  illuminates  the  photocell  and  its 
resistance  diminishes.   This  applies  a  slowly  varying 
positive  voltage  to  the  anode  of  tunnel  diode,  Dl,  through 
an  input  sensitivity  adjustment,  Rl .   When  the  peak  point 
of  the  tunnel  diode  is  reached  it  rapidly  changes  state, 
and  when  it  does,  transistor  Tl  is  driven  into  saturation 
with  a  rise  time  O.lys.   This  pulse  is  now  inverted  by 
transistor  T2  and  both  of  them  are  capacitor  coupled  to 
switch  #1  which  is  used  to  select  either  the  positive 
or  negative  10V  pulse  for  the  address  reset  output  to  the 
multichannel  scaler.   The  positive  pulse  is  also  used  to 
trigger  the  bistable  for  the  constant  velocity  generator, 
and  it  is  brought  out  for  use  on  the  logic  board  on  pin  1 . 
The  constant  velocity  generator  consists  of  a  bistable  (T3, 
T4)  and  a  DC  level  shifter  circuit  (T5,  R16-R19).   T5  inverts 
the  output  of  one  side  of  the  bistable  and  shifts  its  level 
to  +12V.   This  signal  is  now  summed  with  the  complementary 
output  of  the  bistable  through  potentiometer  R19,  and  by 
adjusting  R19  the  voltage  can  be  made  symmetrical  about  any 
DC  level  contained  within  ±12V  with  a  constant  12V  range. 
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c .   Logic  Circuit 

The  logic  circuit  was  designed  and  built  to  inter- 
face the  data  signals  to  the  multichannel  scaler,  to  provide 
a  safety  interlock  circuit,  and  to  gate  the  data  signals  for 
an  integral  number  of  sweeps.   The  circuit  is  shown  in  figure 
61. 

The  data  input  signal  (either  a  pulse  from  the 
detector  system  or  a  pulse  from  the  fringe  discriminator) 
enters  the  board  on  pin  5.   It  is  coupled  to  a  monostable 
multivibrator  by  capacitor  C3  and  coupled  into  the  data  gate 
by  capacitor  CI.   The  monostable  shapes  the  pulse  so  that 
they  are  suitable  for  controlling  the  delayed  coincidence 
gate  on  our  multichannel  pulse-height  analyzer.  This  feature 
allows  one  to  accurately  set  the  single  channel  analyzer 
on  the  peak  in  the  spectrum  of  interest  for  MSssbauer  studies . 
The  data  gate  generates  pulses  to  the  time  mode  data  input 
of  the  multichannel  analyzer  only  under  the  following  con- 
ditions: (1)  the  drive  is  running,  (2)  the  detector  gate 
switch  is  held  in  the  enable  position  until  the  control 
bistable  (T7,  T8)  is  set.   The  control  bistable  is  set  by 
the  first  photocell  trigger  pulse  that  occurs  while  the 
detector  gate  switch  is  held  in  the  enable  position.  When 
the  detector  gate  switch  is  held  in  the  disable  position, 
the  control  bistable  will  be  reset  by  the  next  photocell 
trigger  pulse,  i.e.  the  data  input  signals  are  only  counted 
for  an  integral  number  of  cycles  of  the  Mossbauer  drive. 

The  interlock  circuit  protects  the  drive  against 
photocell  lamp  burnout  and  it  also  interlocks  the  spectrometer 
to  the  analyzer  so  that  the  analyzer  must  be  in  the  time 
mode  store  program  before  the  power  amplifier  of  the  drive 
will  operate.   The  interlock  circuit  consists  of  an  "and" 
gate  and  a  relay.   The  relay  controls  the  DC  bias  to  the 
power  amplifier  so  that  all  of  the  conditions  for  an  "and" 
condition  must  be  present  before  the  relay  will  be  activated. 
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The  conditions  for  an  "and"  are:  (1)  the  analyzer  must  be 
in  the  time  mode  store  program,  (2)  the  photocell  lamp  must 
not  be  burned  out  and  (3)  the  power  switch  on  the  front  panel 
must  be  on.   The  interlock  circuit  also  must  be  in  a  "go" 
position  for  the  gate  to  be  enabled,  so  that  if  the  lamp 
burns  out  during  a  spectrum  accumulation  the  storage  of  bad 
data  will  not  occur.   The  condition  of  the  gate  is  indicated 
by  a  front  panel  lamp  which  is  on  when  all  of  the  data  gate 
conditions  are  satisfied. 

d .   Physical  Description 

(1 )  Layouts 

The  physical  layout  of  components  on  the  chasis 
is  illustrated  in  Figure  62  a,b,c.   Figure  62a  shows  a  layout 
of  the  front  panel.   Figure  62b  shows  the  layout  of  the  chasis 
and  Figure  62c  shows  the  layout  of  the  back  rear  panel. 

The  next  3  figures,  63-5 ,    show  the  layout  of  the 
printed  circuit  cards. 

( 2 )  Wiring  List  See  Appendix  II 

The  wiring  list  is  one  which  is  redundant , 
but  because  of  its  redundancy  it  is  a  self  check  of  almost 
all  of  the  wiring.   The  unit  is  built  with  a  one  point  ground 
system  to  avoid  ground  loop  problems. 

(3)  A  parts  list  can  be  made  available 
upon  request . 

3 •   Frequency  Ratio  Detector 

a.   Introduction 

It  was  decided  to  use  the  information  generated  by 
the  interferometer  in  a  feedback  loop  to  control  the  drive's 
velocity.   Because  the  fringe  frequency  contains  the  velocity 
information,  it  was  necessary  to  compare  the  fringe  frequency 
with  a  known  frequency  and  to  provide  a  signal  which  is 
suitable  for  use  in  a  negative  feedback  loop.   The  frequency 
ratio  detector  was  designed  and  built  to  fulfill  this  need. 
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b .   Description 

The  logic  of  this  device  is  described  elsewhere  l39l , 
and  it  is  only  necessary  to  note  some  of  the  modifications 
which  were  required  for  our  particular  problem.   The  fringe 
detection  system  does  not  determine  the  direction  of  the 
velocity  vector,  and  therefore  the  only  information  obtainable 
from  the  fringe  frequency  (Pf)  is  the  magnitude  of  the  velocity 
or  Ff  a  |v| .  Because  the  sign  of  the  velocity  changes,  it  is 
necessary  to  change  the  sign  of  the  correction  signal,  and 
this  is  accomplished  by  interchanging  the  two  inputs  to  the 
frequency  ratio  detector.   The  inputs  are  interchanged  at  the 
time  of  transition  of  the  "100s"  bistable  in  the  multichannel 
analyzer,  and  the  controlling  signal  enters  the  frequency 
ratio  detector  through  the  velocity  sense  input . 

The  output  of  the  frequency  ratio  detector  is  integrated, 
and  then  summed  with  the  analog  of  the  address  scaler  to  pro- 
duce the  negative  feedback  control  of  the  spectrometer's 
velocity . 

k .   Control  Logic  for  Frequency  Synthesizer 

a .  Introduction 

Because  the  frequency  ratio  detector  needs  a  source  of 
frequency  with  which  to  compare  the  fringe  frequency,  it 
was  necessary  to  design  and  build  a  digitally  programmable 
frequency  source.   This  was  done  by  interfacing  a  remotely 
programmable  frequency  synthesizer  to  the  digital  address 

scaler . 

b .  Interface  Description 

In  Figure  66  a  block  diagram  of  the  control  logic  is 
shown.   The  binary  coded  decimal  (BCD)  address  scaler  infor- 
mation is  decoded  into  decimal  by  the  digital  decoder.   The 
decimal  information  is  then  put  into  a  circuit  which  either 
transfers  the  numerical  information  unchanged  or  takes  its  9Ts 
complement.   For  example  if  the  number  were  28  it  would  be 
transferred  as  a  71  in  the  9's  complement  mode,  or  as  a  28 
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in  the  unchanged  mode.   It  is  necessary  to  complement  the 
channel  address  number  in  channels  0-99  so  that  the  synthe- 
sizer's frequency  is  a  decreasing  function  vs  time,  as  is 
the  fringe  frequency.   The  address  number  should  be  passed 
unchanged  in  channels  100-199  so  that  the  synthesizer's 
frequency  is  an  increasing  function  vs  time,  as  is  the 
fringe  frequency.   The  circuit  which  complements  the  address 
information  is  controlled  by  the  analyzer's  100  bistable. 
The  numerical  information  is  then  transferred  to  the  synthe- 
sizer's remotely  programmable  digital  units  by  the  synthesizer 
driver . 

Performance  of  this  spectrometer  will  be  described  in 
detail  in  a  manuscript  to  be  submitted  for  publication  and 
in  next  years  progress  report . 

(F.  C.  Ruegg  and  J.  J.  Spijkerman) 
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5.   RADIOISOTOPE  TECHNIQUES 

A .   Introduction 

Radioisotope  techniques  have  been  used  in  many  analytical 
techniques  to  provide  data  on  the  adequacy  of  separation  pro- 
cedures and  to  generate  quantitative  results .   Quite  often 
the  results  provide  valuable  information  which  is  difficult 
to  obtain  in  other  ways .  A  considerable  enhancement  to  the 
value  can  be  realized  if  the  precision  of  such  tracer  exper- 
iments could  be  improved.   It  is  in  this  area  that  our  atten- 
tion is  being  focussed. 

Reproducibility  of  the  radiation  detection  process  is 
being  evaluated  along  with  the  chemical  techniques  to  deter- 
mine which  variables  are  limiting.   In  particular,  for 
macro-constituent  elemental  quantitative  analysis  it  is 
important  to  utilize  reproducible  chemical  reactions  in  the 
separation  process.   The  substoichiometric  approach  may  be  a 
good  one  for  improving  this  precision.   This  may  be  the  case 
particularly  when  interfering  chemical  reactions  are  likely. 

In  substoichiometric  radioisotope  dilution  analysis 
(SRDA)  (4l-2),  the  need  to  isolate  the  element  of  interest 
quantitatively  (a  requirement  of  most  of  the  other  analyti- 
cal techniques)  is  eliminated  and  the  determination  of  the 
specific  activities,  before  and  after  dilution  by  the 
analyzed  element,  can  be  replaced  by  activities  in  the 
relation : 


X  =  Y 


'1 
■  -  1 


LA2 


Where  X  =  the  amount  of  element  to  be  determined;  Y  =  the 
known  amount  of  element  in  the  standard  solution;  A-.  =  the 
radioactivity  of  the  standard  solution  and  A„  =  the  activity 
of  the  unknown  solution.   The  activities  A.  and  A„  are 
usually  obtained  as  counting  rates  with  a  suitable  elec- 
tronic counter,  thus  enabling  one  to  associate  a  minimum 

relative  error  with  the  counting  rate  which  would  be 
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affected  only  by  counting  statistics  if  experimental  con- 
ditions are  kept  constant  (i.e.,  constant  geometry  maintained 
and  repeated  measurements  made  [43]).  In  most  cases,  the 
error  associated  with  Y  is  known  and  very  small,  therefore, 
one  should  be  able  to  do  chemical  analysis  (determine  X) 
very  accurately  and  precisely  with  the  substoichiometric 
radioisotope  dilution  technique. 

The  first  part  of  this  report  describes  the  counting 
system  used  for  determining  the  activities  A,  and  A? 
mentioned  earlier  and  the  precautions  taken  to  yield  the 
minimum  error  obtainable.   The  latter  part  of  this  report 
describes  applications  of  the  SRDA  technique  to  several 
elements . 
B .   Error  Analysis  of  Counting  System 

1 .  Instrumentation 

The  instrument  used  for  our  work  is  very  simple  and 
economical,  which  is  a  great  advantage  to  low  budget  insti- 
tutions interested  in  high  precision  quantitative  analysis. 
It  consists  of  a  2in.  x  2in.  Nal(TI)  well-type  scintillation 
detector,  mounted  very  rigidly  so  there  would  be  no  vibration, 
enclosed  in  a  four  inch  thick  shield  built  with  lead  bricks 
and  a  single  channel  pulse  height  analyzer.  The  use  of  this 
quality  of  shielding  is  not  critical  at  the  count  rates  used. 
The  system  is  used  for  other  counting  purposes  in  addition 
to  this  one.   A  plastic  liner  which  was  made  to  fit  firmly 
into  the  detector  serves  a  dual  purpose,  to  hold  the  counting 
vessels  (vials  which  are  screened  for  uniformity)  in  a  con- 
stant geometry  and  to  prevent  any  contamination  of  the  detec- 
tor.  Photographs  of  the  detector  and  its  accessories  are 
shown  in  Figures  67  and  68.   The  cost  of  this  system  is  less 
than  thirty-five  hundred  dollars . 

2 .  Calibration 

Counting  the  same  volume  of  sample  solution  each  time 
allows  one  to  maintain  the  constant  geometry  with  the  system 
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Figure  67.   Photograph  of  radiation  detector  assembly. 


previously  described,  so  that  the  error  in  the  counting  rates 
(A-.  and  A?)  will  be  due  mainly  to  counting  statistics. 
Cobalt-60  was  chosen  as  the  standard  isotope  to  calibrate  and 
determine  the  performance  of  this  counting  system.   The  con- 
trols are  adjusted  on  the  analyzer  to  exclude  bremsstrahlung 
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Figure  68.   Photograph  of  electronic  apparatus  medium 
count-rate  capability. 


so  that  only  the  photopeak  of  interest  is  counted.   This 
reduces  the  sensitivity  of  the  system  to  source-position 
errors.   Counting  an  aliquot  of   Co  repeatedly  and  collecting 
at  least  4  x  10   counts,  resulted  in  a  reproducibility  of 
0.04%  relative  standard  deviation  of  a  single  measurement. 
Further  study  is  underway,  utilizing  high  speed  electronic 
equipment  (scaler,  timer,  etc.)  which  may  improve  our  pre- 
cision by  a  factor  of  ten. 
C .   Application  to  Standard  Reference  Materials  (SRM) 

Our  research  efforts  have  been  directed  to  those  elements 
in  which  this  technique  would  be  superior  to  other  existing 
quantitative  analytical  techniques  (e.g.  spectrophotometry, 
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gravimetry)  and  thus  play  a  vital  part  in  the  SRM  certifi- 
cation program  at  the  National  Bureau  of  Standards.   At  the 
same  time  a  high  degree  of  flexibility  is  maintained  so  that 
the  method  developed  could  be  adapted  easily  for  both  routine 
analyses  and  special  cases. 

Before  an  element  is  considered  a  candidate  for  analysis 
by  the  SRDA  technique,  it  must  have  a  long-lived  radioisotope 
obtainable  with  a  high  specific  activity.   Table  7  lists  those 
elements  where  use  of  the  SRDA  technique  would  be  advantageous 
not  only  to  the  Bureau's  SRM  program,  but  to  industrial  labor- 
atories as  well,  and  a  significant  contribution  to  the  state 
of  the  art  is  accomplished  simultaneously  with  the  development 
of  these  methods. 

Table  7.   Elements  that  can  be  analyzed  by  the  substoichio- 
metric  radioisotope  dilution  analysis  technique 
in  the  Standard  Reference  Materials  Program. 

Antimony  Niobium 

Cadmium  Silver 

Chromium  Tin 

Cobalt  Titanium 

Iron  Tungsten 

Mercury  Zinc 


Previous  reports  [3-^1]  described  the  use  of  the  SRDA 
technique  for  trace  elemental  analysis  (Co,  Cd,  Ag,  etc.) 
in  several  Standard  Reference  Materials.  In  this  report  we 
will  describe  the  SRDA  methods  developed  for  tungsten  and 
iron  and  discuss  the  results  obtained  when  SRM  480  and  SRM 
1159  and  1160  were  analyzed  for  these  elements  using  this 
technique . 
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1 .   Determination  of  Tungsten  in  Molybdenum  Containing 

Materials  Using  Toluene-3 ,4-Dithiol 
There  are  several  reagents  which  produce  colored  complexes 
with  tungsten  (44)  thus  providing  the  basis  for  numerous 
spectrophotometric  methods  for  the  determination  of  tungsten 
in  various  materials.   However,  these  reagents  also  produce 
colored  complexes  with  molybdenum,  and  since  most  tungsten 
containing  materials  also  contain  molybdenum,  an  interference 
problem  exists  with  these  methods.   Most  colorimetric  methods 
which  result  in  the  formation  of  stoichiometric  complexes 
can  be  adapted  to  radioisotope  dilution  techniques. 

Toluene-3 ,4-dithiol  forms  colored  complexes  with  several 
elements  [45]  including  tungsten  and  molybdenum  in  various 
materials.   Several  workers  [46-50]  have  developed  colorimetric 
methods  for  determining  tungsten  and  molybdenum  but  the  success 
of  these  methods  is  very  dependent  upon  the  control  of  the 
problems  (quantitative  separation  of  each  element,  interfer- 
ences, color  stability,  etc.)  inherent  in  photometric  methods. 
Thus  arising  from  the  need  of  a  method  in  our  laboratory  for 
tungsten  determinations  in  the  presence  of  molybdenum,  we 
developed  the  following  radioisotope  dilution  procedure  for 
tungsten  using  "dithiol" . 
a .   Experimental 

The  analyzed  samples  were  special  cuts  from  NBS  W-Mo 
Alloy  (SRM-480).   Each  sample  (^5  mg)  was  etched  in  a  freshly 
prepared  mixture  of  10  ml  concentrated  HNO    10  ml  cone. 
H2SOi[  and  10  ml   H20  for  15  minutes.   The  samples  were  rinsed, 
dried  and  weighed.   The  weight  loss  from  this  etching  step 
is  very  small  (<  .5%)  .      Each  sample  was  dissolved  in  approx- 
imately 10  drops  HF  plus  3  drops  HNO    heated  until  NO   was 
driven  off,  diluted  to  5  ml  in  a  volumetric  flask  with  1  M 
HC1  and  then  transferred  immediately  to  a  polyethylene  bottle. 
Aliquots  of  these  samples  were  taken  for  analysis  .   Prior  to 
tungsten  analysis,  a  molybdenum  separation  was  carried  out 
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as  follows:  8.5  ml  HC1  (sp.  gr .  1.06),  a  known  amount  of  W 

-i  O  -i 

carrier  (Y)     W  and  4  drops  of  10%  NH20H.HC1  were  added  to 
the  sample  and  mixed.   Then  0.15   ml  of  0.5%  dithiol  solution 
was  added.   After  5  min  the  Mo-dithiol  complex  was  extracted 
into  CCIk  (four  times  using  5  ml  portions).   The  aqueous 
phase  (~10  ml)  was  evaporated  very  carefully  to  2  ml  and 
25  ml  hot  titanous  sulfate  and  20  ml  hot  cone.  HC1  were 
added.  To  the  mixed  solution  0.20  ml  of  the  dithiol  solution 
was  added,  mixed  and  allowed  to  stand  for  90  minutes,  after 
which  the  W-dithiol  complex  was  extracted  into  10  ml  CHCl^, 
and  0.50  ml  aliquots  of  both  phases  were  taken  for  counting. 
A  control  (i.e.,  all  ingredients  added  except  the  sample) 
is  subjected  to  the  same  procedure,  providing  A-.  . 

If  no  molybdenum  is  present,  the  Mo-dithiol  complex 
formation  can  be  omitted,  and  the  remaining  part  of  the  pro- 
cedure would  be  carried  out  as  described.   The  amount  of 
tungsten  present  is  calculated  using  equation  32. 

b .   Results  and  Discussion 

The  results  obtained  by  this  radioisotope  dilution 
method  are  reported  in  a  note  to  be  published  in  the  journal, 
Analytical  Chemistry.   The  following  is  an  abstract  of 
this  article : 

RADIOISOTOPE  DILUTION  ANALYSIS  OF  TUNGSTEN  AS  A  MAJOR 
CONSTITUENT  IN  MOLYBDENUM  CONTAINING  MATERIALS 
USING  TOLUENE-3,4-DITHIOL 
ABSTRACT 
A  radioisotope  dilution  procedure  was  developed  for 
tungsten  determination  in  molybdenum  containing  materials 
using  a  substoichiometric  amount  of  the  complexing  reagent, 
toluene-3, 4-di thiol.   Using  this  approach,  an  equal  amount 
of  tungsten  is  isolated  before  and  after  dilution  of  the 
radioisotope,  thus  avoiding  serious  interferences  by  elim- 
inating the  need  for  quantitative  separation.   This  tech- 
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nique  was  used  for  tungsten  determination  in  NBS  tungsten- 
molybdenum,  where  the  average  concentration  found  was 
79 . 31%  with  a  pooled  estimate  of  the  standard  deviation 
(for  a  single  determination)  of  0.43%.   The  precision 
obtained  is  very  good  for  this  technique . 

2.   Determination  of  Iron  in  SRM  1159  and  1160 
Iron  is  one  of  the  elements  which  is  very  often  certi- 
fied at  varying  concentrations  in  (SRM)  steel  samples  at 
the  Bureau  of  Standards.   A  radioisotope  dilution  procedure 
was  developed  for  iron  which  utilizes  solvent  extraction  [51]. 
This  procedure  offers  an  attractive  alternative  to  existing 
chemical  methods  for  the  determination  of  iron  because  it 
provides  a  rapid  analysis  with  very  good  precision  and  can 
be  adapted  to  iron  concentrations  at  all  levels .   Iron  was 
determined  in  SRM  1159  where  it  is  the  major  constituent  and 
in  SRM  ll60  where  it  is  a  minor  constituent  in  a  nickel  matrix 

a.  Experimental 
— ± — 

Samples  of  these  SRM  steels  weighing  one  gram 
were  dissolved  in  concentrated  HC1-HN0   mixtures  (minimum 
amount)  and  diluted  to  100  ml  in  volumetric  flasks  with  H  0. 
Aliquots  (1  ml  -  3  ml)  of  the  samples  were  placed  in  separa- 
tory  funnels  along  with  a  known  amount  of  Pe+3 (Y) ,59Fe , 
3  ml  of  30%  H202  and  mixed  well.   0.75ml  of  2%    cupferron  was 
added  to  each,  mixed,  and  the  Fe-cupferron  complex  extracted 
with  10  ml  CHCl^.   The  volume  of  the  aqueous  phase  is  adjusted 
to  10  ml  with  water  before  extraction.   A  control  (A  )  was 
subjected  to  the  same  procedure.   Aliquots  (0.5  ml)  of  both 
phases  were  taken  for  counting  and  the  amount  of  iron  present 
was  calculated  using  equation  32. 

b .  Results  and  Discussion 

The  results  of  the  iron  analysis  are  shown  in  Table 
8.   The  precision  obtained  with  this  technique  is  excellent. 
Precautions  must  be  taken  to  insure  that  all  of  the  iron 
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present  is  in  the  ferric  state.   For  this  reason  hydrogen 
peroxide  is  added  in  the  procedure . 

Iron  is  an  element  that  is  of  biological  interest. 
Thus,  this  technique  can  be  used  for  iron  determinations  in 
varied  matrices  with  rapidity  and  very  low  cost. 

Table  8.   Analysis  of  Iron  in  Standard  Reference  Materials 

SRM  %   Fe  Average 

1159  51.21  51.19i0.0081 

51.19 

51.19 
51.18 
51.20 
51.19 
51.17 
51.19 

1160  14.51  I4.50t0.0081 

14.49 

14.50 
14.51 
14.51 
14.50 
14.50 


t  s 
1   Errors  are  given  as : at  the  95%  confidence  level  [52] 
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Appendix  I.   Evaluation  of  Integrals 

Evaluation  of  the  Integrals  <n, I ,m| 0 | n' , I " ,m'> 


0  = 


8  - 


A,    =    /°°dy    yJ"    exp[-    ^(^  +   §C)p] 


J    i 


2  n    n 


B.  =  /  dv  vJ"  exp[-  i(^i  +  B|l)v] 
j    -i  ^    2Nn    n 


[2n+l]/2 


C   =  (R?) 


[2n'+l]/2 


D   =  (R^) 


<1,0,0 

o| 

/N 

<1,0,0 

0| 

/s 

<1,0,0 

o| 

^ 

<1,0,0 

o| 

/N 

<1,0,0 

o| 

•s. 

<1,0,0 

o| 

1,0, 0>   = 


\      I  CD^AiBo+AoBi] 


2>°>0>   =   M   R  CD[A2B0+2AlBl+A0B2] 


350,0>   =  32^3Q  |  CD[A3B0  +  3A2B1+3A1B2  +  A0B3] 


2,1,0>   =   j^  a  CD[A2B1+A1(B2+B0)+A0B1] 


3,1,0>      =    32Wi0RCD^A3Bi+A2(2B2+   BQ  HA-^B^B^+AqB^ 


3'2'0>   =  SW^   R  CD[A3(3B2-B0)  +  3A2(B3+B1) 

+3A1(B2+B0)+A0(-B3+3B0)]. 
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<2,0,0 | 0 | 2,0,0>  = 


lh        R  CD[A3B0  +  3A2B1+3AlB2  +  A0B3] 


1 Q 


L2,0,0|O|3,0,0>    =    155^2t/5   g   CD[A4B0+4A3B1+6A2B2+4A1B3+A0B4] 

:2,0,0|0|2,1,0>    =         y|g-        |   CD[A3B1+A2(2B2+B0)+A1(B    +2B1) 

+a0b2]. 


q 


2,0,0|O|3,l,0>    =    5ie^15   g   CD[Ai|B1+A3(3B2+BQ)+3A2(B3+B1) 

+A1(B4+3B2)+AQB3]. 

<250s0|0|3,2,0>    =      31^Qi|      |  CD[A4(3B2-B0)+A3(6B3+2B1) 

+A2(3B4+6B2+3Bq)+A1(2B3+6B1) 
+A0(-B4+3B2)]. 


q 


<3,o,o|o|3,o,o>  =  24o78ia  r  cdu^o+sa^+ioa^+ioa^ 

+  5A1B4+AQB5].. 

<3,0,0|6|2,1}0>    =    51q^15':§  CDLA^B^    A3  (  3B2+BQ  )+A2  (  3B3  +  3B1) 

+A1(B1|+3B2)+A0B3]. 

<3,0,0|o|3,l,0>    =    2i|Q^8l^    |  CD[A5B1+A1|(4B2+B0)+A3(6B3+4B1) 

+a2(4b4+6b2)+a1(b5+^ib3  )+a0b4] 
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<3,050|0|3,2,0>    =    it8^8l2    |   CD[A5(3B2-B0)+A4(9B3+B1) 


+  A3(9Blj+8B2+3B0)+A2(3B5+8B3+9B1) 
+A1(B1|+9B2)+A0(-B5+3B3)]. 


1      Q 


<2,1,0|0|2,1,0>    =      2%  r    CD[A3B2+A2(B3+2B1)+A1(2B2+B0)+A0B1] 


-i_  ^ 


<2,l,0|O|3,l,0>    =    5YM75    R    CD[A4B2+A3(2B3+2B1)+A2(Bi|  +  4B2+B0) 


+A1(2B_+2B1)+A0B2]. 


<2 ,1,0 | 0 | 3,230>  = 


10368/3  R  CD[A4(3B3-Bl)+A3(3B4+6B2-B0) 

+A2(6B-+6B1)+A1(-B1.+6B2+3B0) 


+AQ(-B  +3B1)]. 


<3,1,0.|0|3,1,0>    =      8o.8l^    |   CD[A5B2+A4(3B3+2B1) 


+A3(3Bij  +  6B2+B0)+A2(B5+6B3+3B1) 
+A1(2B2|+3B2)+A()B3]. 


<3,1,0|0|3,2,0>   =    1|8^T|1.   |  CD[A5(3B3-B1)  +  A4(6B4+5B2-B0) 


+A3(3B5+12B3+5B1)+A2(5Blj+12B2+3B0) 


+  A1(-B5+5B3+6B1)+AQ(-B1|+3B2)] 
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<3,2,0|0|3,2,0>  =  192,8l2  |  CD[A5(9B4-6B2+B0) 

+Ai|(9B5+l8B3-7B1) 

+A  (18b^+28B2-6bq) 

+An(-6B[-+28B_  +  l8Bn  ) 
2     D  5  1 

+A1(-7Bi|  +  l8B2  +  9B0) 
+A0(B5-6B3+9B1)]. 

f2,l,±l|0|2,l,±l>  =  5X2  §  CD[(A3-A1)(B0-B2)+(A2-A0)(B1-B  )] 


2,l9±l|0|3,l,±l>    =    103^8v/5   I  CD[(A4-A2)(BQ-B2) 

+  2(A3-A1)(B1-B3)  +  (A2-A0)(B2-Bi|)] 


1    q 


2}13±1|0|3S23±1>  =  Iq|m  ^CD[(Ai|-A0)(B1-B3)  +  (A3-A1)(B0-B4)] 


1 q 


3,1,±1|0|3,1,±1>  =  1g0:8l2gCD[(A5-A3)(B0-B2)  +  3(Alj-A2)(B1-B3 


+  3(A3-A1)(B2-B]|)  +  (A2-A0)(B3-B5)]. 


<3,l,±l|0|3,2J±l>  =  l6^8l.  |CD[(A5-A1)(B1-B3)  +  (Ai|-A2)(B0-B1|) 


<3}2}±1|0|352,±1>  = 


<3i25±2|0|3,2,±2>  = 


+(A3-A1)(B1-B5)+(A4-AQ)(B2-B4)]. 

32.3X2  I  CD[(A5-A1)(B2-Bij)  +  (A4-A2)(B1-B5) 
+  (A3-A1)(B0-B1|)  +  (A4-A0)(B1-B3)] 

Ig^ep-^  CD[(A5-2A3+A1)(B0-2B2+B4) 

+(A4-2A2+AQ)(B1-2^  +B5)]. 
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Appendix  II 


Pin  # 


2,  14,  18 

3 
5 


6,  24 


12 
13 
15 
16 
17 
21 
22 
23 


Wiring  List  for  MSssbauer  Spectrometer 
Board  #1 
Signal 


Detection  pulses 
to  analyzer 

-12V 


Ground 
Data  Signal 

+  12V 


Photocell  Lamp 
Burnout  interlock 


To 

S.C.  #2  R 
Pin  #3 

Power  Plug 
Pin  #17 

Chasis  Grnd .  Buss 

Center  Pole 
Switch  §3   F 

Power  Plug 
Pin  #16 

S.C.  #1 

Pin  #2 


To  delay  coincidence  S.C.  #2 
input  of  analyzer     Pin  #2  R 


Relay  controlled 
-12V  =  -12R 

Detector  gate 
Lamp  Driver 

Relay  Controlled 
+12V  =  +12R 


Bd.  #3 
Pin  #3 

Lamp  F 


Bd  #3 

Pin  #17 


Analyzer  Time  Store   S.C.  #2   R 


Interlock  input 

Control  Bistable 
gate  line  (enable) 

Position  Photocell 
Trigger  Pulse 

Control  Bistable 
gate  line  (inhibit) 


Rn  #1 

Switch  #4  P 
Bottom  Pole 

Bd  #2 
Pin  #1 

Switch  #4  F 
Top  Pole 


128 


Mossbauer  Photocell  Trigger 
Board  #2 


Pin  #         Signal 


Positive  photocell 
Trigger  pulse 

Address  reset 


3  +12V 

4  Photocell 
Bias  +12V 

5  Photocell 
Signal 

7  Ground 

11  -12V 


19  Positive 

Square  Wave 

21  Negative 

Square  Wave 


Wire  Or 
Coax 

To 

W 

Bd  #1 
Pin  #22 

w 

Switch  #2  F 
Center  Pole 

w 

Bd  #1 
Pin  #6 

w 

S.C.  #1  R 
Pin  #8 

w 

S.C.  #1  R 

Pin  #9 

w 

Ground  Buss 

w 

Bd  #1 
Pin  #14 

c 

Pot  #3  F 
C.  W.  Pole 

c 

Pot  #3  F 
C.C.W.  Pole 
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Mossbauer  Amplifier 
Board  #3 


Wire  or 
Pin  §  Signal  Coax        To 

■ 
1  +12V  W         Power  Plug 

Pin  #16 

2,6,9,         Ground  W         Ground  Bus 

16,  22 

3  -12V  W         Power  Plug 

Pin  #17 

4  Input  to  C 
Spare  Op .  amp 

5  Output  of  C         BNC  F 
Spare  op.  amp 

7  L.V.  Feedback  C         S.C.  #1  R 

Signal  Center  Pole 

12  C         Base  of 

STC  1080 

13  -12V  R  W         emitter  of 

STC  1080 
and  Bd  #1, 
Pin  #13 

15  Power  amp  C         Collectors 

output  of  STC  1080 

and  2N3171  & 
BNC  R 
Power  out 

17  +12V  R  W         emitter  of 

2N3171  and 
Bd  #1,  Pin  #16 

19  C         Base  of 

2N3171 

20  Symmetry  C         Pot  #2 
Control  Center  Pole 
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F 


Pin   § 


23 
24 


Wire  or 

Signal 

Coax 

To 

Address  Analog 

C 

S.C 

.  #2  R 

Input 

Pin 

#6 

Input  Signal 

C 

Pot 

#1  P  Pin  #3 

Attenuator 

BNC 

P  Input 

Square  wave 

C 

Pot 

#3  F 
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Signal  Connector  #1  (S.C.#1)  Mounted  on  Rear 

Wire  or 

Pin  #  Signal  Coax  To 

2  Photocell  lamp           W  Bd  #1 
burnout  interlock  Pin  #7 

3  Photocell  lamp           W  Switch  #1  F 
Power  Bottom  Pole 

8  Photocell                W  Bd  #2 
Bias  +12V  Pin  #5 

9  Photocell                W  Bd  #2 
Signal  Pin  #5 

11  Ground  W  Ground  Buss 

13  L.V.  Feedback         Center  Bd  #3 
Signal                Conductor  Pin  #7 

of  coax 

14  L.V.  Feedback         Shield  of  Ground 
Signal               above  coax   at  Bd  #3 
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Signal  Connector  Cable  #1  to  Drive  Motor 

All  Coax 

Pin  #  Signal  To 

2  Lamp  interlock  Lamp 

3  Lamp  power  Lamp 

8  Photocell  Bias  Photocell 

9  Photocell  signal  Photocell 

11  Ground  Ground 

13  Linear  Velocity  LV 

Transducer  Feedback 
Signal 

14  Linear  Velocity  LV 
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S.C.#2  Mounted  on  Rear 


Wire  or 

Pin  #                       Signal  Cable        To 

1  Analyzer  Store  W         Bd  #1 
interlock  Pin  #17 

2  To  delay  coincidence  W         Bd  #1 
input  of  analyzer  Pin  #12 

3  Time  Mode  Detector  W         Bd  #1 
input  of  analyzer  Pin  #1 

4  To  output  of  Single  W         Switch  #3 
Channel  Analyzer  Bottom  PoleP 

6           Analog  of  Address  C         Bd  #3 

Scaler  Pin  #21 

8           Address  Reset  W         Bottom  Pole 

Sw  #2 

10  To  fringe  pulse  Switch  #3  F 
train  Top  Pole 

11  Ground  W         Ground  Buss 
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Signal  Connector  Cable  #2 

All  coax  with  BNC  coaxial  connector  on  other  end  of  each 
connection. 


Pin  #        Signal 


To 


Connector 


Analyzer  store 
Interlock 


Multichannel 
Analyzer 


P-2  Pin  J 

External  Live  Timer 


Delay  Coincidence  Multichannel 
Pulses  Analyzer 


Delay  Coincidence 


3 

Data  Signal 

Multichannel 
Analyzer 

Time  Mode  Detector 
Input 

k 

Output  of  Singal 
Channel  Analyzer 

Single 

Channel  Analyzer 

Output 

6 

Analog  of  Address 
Scaler 

Multichannel 
Analyzer 

Ramp  Out 

8 

Address  Reset 

Multichannel 
Analyzer 

Address  Reset 
Input 

0 

Fringe  count 
Pulse  train 

Frequency 
Ratio  Detector 

Fringe  Pulses 
Out 

11 


Ground 


Shield  of 

All  Coax  Cables 
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Pot  #1 


Pot  #2 


Pot  #3 


Front  Panel 
Velocity  Amplitude  Control 
use  coax  20  K 

Pin  #3 


10  turn 


Pin  #2 


Pin  #1 


BcT  #3 
Pin  8 


/TnBNC® 
v_y  Input 


vBd  #3 


Pin  2  3 


Constant  Acceleration  Symmetry  Control 
use  coax  on  wiper     100  K    10  turn 


Bd  #3 


Pin  #3 

+  12V 

XPin  #2                                      J 

^         — 

r Pin  #1 

r 

-12V 

Pin  20 


Constant  Velocity  Symmetry  Control 
Use  coax  10  K    one  turn 


CW 


Bd  #2 
Pin  19 


CCW 


vBd  #3 
Pin  24 


.Bd  #2 


Pin  21 
Switch  #1   Power  Switch 


on 

/ 

S.C.  #1 
Tin  #3 

-12V 

off 

Switch  #2 

Address  Reset 

S.C.  #2 
,               Pin  #8 

/ 

Bd  #2 
Pin  #2 


136 


Switch  #3   Mossbauer  or  Calibrate 
Mossbauer 

, JS.C.  #2 

Pin  #4 


Bd  #1 
Tin  #5 


, S.C.  #2 

Pin  #10 
Calibrate 

Switch  #4   Control  Bistable  Gate 


Store Bd  #1 

Tin  #21 


-ll' 


Bd  #1 
Pin  #23 


Inhibit 
Lamp  Indicator  for  Detector  Gate 


■12V fe§ *Bd  #1 

Tin  #15 


BNC  Output  of  Spare  Op.  Amp. 
Use  Coax 


©■ 


^Bd  #3 
Pin  #5 
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